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When you can measure what you are speaking about, and express it in 
numbers, you know something about it; but when you cannot measure 
it, when you cannot express it in numbers, your knowledge is of a meager 
and unsatisfactory kind: it may be the beginning of knowledge, but you 
have scarcely, in your thoughts, advanced to the stag^ of science. 

- Sir William Thomson, Lord Kelvin of Largs (1824-1907) 
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ABSTRACT 

A hydro6t&tic primitive equation model initialized in a highly baroclinicaUy unsta¬ 
ble state was used to simulate maritime cyclogenesis and frontogenesis. In order 
* to identify boundary layer physical processes important in maritime frontogenesis, 
several different simulations woe performed. An adiabatic and inviscid simulation 
provided the control for these experinoents. The two different boundary layer pa- 
rameterizations used were a K-theory parameterization and a second-order closure 
scheme. Results indicated that strong warm and cold frtmts formed in the adiabatic 
and inviscid case but that the near-surface wind speed and vertical motion fields 
were unrealistic. In the K-theory simulation, the results were somewhat more realis¬ 
tic but convergence and vorticity were weaker. Results from the second-order closure 
simulation demonstrated that turbulent mixing of momentum was most important 
in producing the frontogenetic (and frontolytic) effects of the transverse secondary 
circulation. 
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I. INTRODUCTION 


Much of the rich variety of sensible weather experienced by human-kind is 
associated with fronts. In many instances, severe weather also is rdated to ascent 
in frontal zones. It is little wonder, therdbre, that fronts have ben a subject of 
investigation since the inception of the polar front theory during the First World 
War. The structure of mature frontal zones was fairly well understood at that time 
due to careful observations and the brilliant insights of Bjerknes and his co-workers. 
The mechanisms responsible for the genesis of fronts, however, were not addressed 
in the frontal/cyclone conceptual model of the Bergn School. Indeed, the theory 
held that the polar front was a semi-permannt, quasi-steady phnomnon whidi 
varied in intensity with occurrences of cyclogenesis along it; fronts were thought to 
be the cause of, rather than a consequence of, cyclogenesis. This concept dianged 
vrith advances in instability theory dturing the 1940’s. Much of the progress in the 
understanding of frontogenesis was mzule using numerical investigations which relied 
on technology not available until the mid 1960*s (e.g., Williams, 1967), although the 
introduction of the frontogenetic function by Petterssen (1936) as extended by Miller 
(1948) represented a major contribution. Investigation of fronts and frontogenesis 
continues to occupy a prominent place in the scientific literature; Keyser and his 
co-workers recently discussed a generalization of the frontogenetic function (Keyser 
et al., 1988) and, more recently, Shapiro and Keyser proposed a major revision to 
the frontal/cydone conceptual model of the Bergen School (Shapiro and Keyser, 
1990). 
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AltbcHigh oonnderable progren haa ben made, some imporUiil qaestions re¬ 
main unanswered. Thne b some oontroveny over the Umiting scale of fronts and 
the existence of a “brake” on frontogenesb due to a separation of the zone of max¬ 
imum convergence from the zone of maximum rdative vortidty. Issues such as the 
• impact of surface heat fluxes on frontogeneus are unres(dved. In fact, some of the 
most important outstan^ng issues concern the impact of boundary layer processes, 
including (but not limited to) surface fluxes oi both heat and momentum. 

The present study constitutes an atteziq>t to address some of the topics related 
to the impact of boundary layer processes on maritime frontogenesb. The approach 
consists of using two different boundary byer parameterization schemes in a numer¬ 
ical simulation of frontogenesis forced by a growing baroclinic wave. An adiabatic 
and invisdd simulation is also performed as a control. The lower boundary condi¬ 
tions are characteristic of an ocean surface. Various processes are isobted in these 
schemes in an effort to understand thdr influence on frontogenesb. In particular, 
surface fluxes of heat and momentum are removed individually from simubtions 
using each of the parameterization schemes. Comparison of these simulations to 
the full physics simulations and to the adiabatic and invisdd simubtion provides 
insights into the impact of surface fluxes on frontogenesb. The ability of the two 
parameterization schemes to realbtically represent the effect of turbulent mixing on 
frontogenesis in general and on the thermally direct ageostrophic drculation about 
the front in particular is also assessed. 

A review of previous work on frontogenesb b presented in Chapter II. Chapter 
III describes the objectives of this study and detaib of the model and numerical 
experiments are discussed in Chapter IV. The results are described in Chapter V 
and dbcussed in Chapter VI. In Chapter VII, the condusions are presented. 
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11. LITERATURE REVIEW 


A. INTRODUCTION 

The introduction of the polar frcmt theory by the Bergen School ibOowing 
the First World War, whidi marked the beginning of meteondc^ as a quantitative 
physical science, hdd that fronts were the dynamic interfaces between differing air 
masses. The suggestion that fronts frnm in growing nonlinear barodinic waves came 
during the 1940*s due to advances in the theory of barodinic instability attributable 
to Chamey (1947) and Eady (1949). This suggestion was further substantiated by 
Phillips (1956) and was demonstrated condusivdy by Williams (1967). Observa¬ 
tional studies have delineated the important structural features of fronts (Sanders, 
1955; Ogura and Portis, 1982; Carbone, 1982; Shapiro, 1984; Fleagle et al., 1988; 
Shapiro and Keyset, 1990), although some controvosy has arisen over the limiting 
scale of frontal zones (Gall et al., 1987). Several analytical and numerical studies 
of the processes by wh^ frontogenesis occurs have been performed (a brief review 
appears in Gill (1982)). However, some important questions regarding frontogenesis 
remain unanswered. One of the most fundamental of these concerns the impact of 
boimdary layer physics on the evolution of fronts. 

B. FRONTS AND BAROCLINIC WAVES 

Although the polar front theory and subsequent studies during the first half of 
this century lead to a fairly detailed description of fronts, the theoretical basis for 
the formation of fronts lagged far behind. The problem which had to be addressed 
was the creation of small-scale features in a field initially containing only large-scale 
features. The energy bearing barodinic modes have scales on the order of 1000 
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km while finmtal scales axe 100 km or less. Williams (1967) used the hydrostatic 
Bottssinesq equations in a numerical model to demonstrate that a frontal discon¬ 
tinuity could be produced in a finite time from an initial disturbance with a scale 
similar to that of the most unstable Eady wave. An analytical treatment of the 
semi-geostrophic equations by Hoskins and Bretherton (1972) (hereinafter, HB72) 
showed similar results. 

Most studies of frontogenesis have focused on the fronts themselves. Less 
attention has been focused on the role of fronts as a coupling agent between large- 
scale baroclinic waves and small-scale dissipative modes. The spectral characteristics 
of semi-geostrophic fronts and baroclinic waves have been investigated by Andrews 
and Hoskins (1978) and Gall et al. (1979), respectively. These studies suggest 
that the spectrum of motions with zonrd wavraumber >10 may be representative 
of processes associated with frontal formation rather than inertial energy transfer 
(Blumen, 1980). 

C. STRUCTURE OF FRONTAL ZONES 

Several years prior to the introduction of the polar front theory, Margules 
(1906) devdoped a detailed theory of atmospheric discontinuities. His work antic¬ 
ipated, in many respects, the dynamical views that were to be introduced some 
40 years later in that he appreciated the importance of the conversion of avulable 
potential energy to kinetic energy. A historical account of subsequent development 
of the current understanding of atmosphwc fronts is given by Palmen and New¬ 
ton (1969). Detailed descriptions of the structure of intense cold fronts appear in 
Sanders (1955) and Petterssen (1956). Typically, cold fronts lie in confluent regions 
associated with a frontal trough. Convective clouds and precipitation accompany 
the front, often in bands. Warm fronts are accompanied by stratiform clouds. Warm 
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front! have recently become the focus of renewed interest. Field observation pro¬ 
grams associated with GALE and ERICA showed that strong warm fronts play an 
important role in cyclogenesis over the western North Atlantic ocean (Nuss, 1989; 
Chang et al., 1991; Neiman et al., 1991). Cold and warm fronts generally intersect 
the center of an extratropical cyclone which is an area diaracterized by extensive 
heavy precipitation. At later stages of frontogene^/cyclogenesis, the site of the 
intersection is the occluded front formed when the cold front overtakes the warm 
front and forces it aloft. On the scale of synoptic diaits, the front is represented 
by a zeroth order discontinuity in wind and a first order discontinuity in pressure. 
Conventional wisdom holds that turbulence and diffusion lead to frontal scales on 
the order of 10-100 km. Recent studies, however, have indicated that some fronts 
may have scales on the order of 100 m (Shapiro et al., 1985; Carbone, 1982). A 
review of observational investigations of fronts is g^ven by Keyser (1986). 

A major revision to the classical Norw<^an frontal cyclone model was recently 
proposed by Shapiro and Keyser (1990). In the Shapiro and Keyser modd, the warm 
and cold fronts are separated near the center of the low by a '%ontal fracture'’. 
Moreover, the warm front extends to the west of the low center; this is referred to 
as the “bent back” warm front. In the later stages of development, a warm core 
“seclusion” forms near the center of the low. There is no occlusion in this model. 
The Shapiro and Keyser model is discussed at length in the Appendix. 

FVontal zones are smallest and frontal gradients most intense at the surface. 
The gradients within the surface frontal zone may be one or two orders of magni¬ 
tude larger than those observed in the free atmosphere (Sanders, 1955). The surface 
frontal zone is characterized by strong horizontal convergence and cyclonic vorticity 
with weak anticyclonic vorticity behind the front associated with the cold air out¬ 
break (Ogura and Portis, 1982). This zone has small vertical extent 100 mb). 
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The front beocmm more diftue eloft and extends thiooi^ the depth ot the tiopo' 
sphere. IWhulent fluxes ci sensible heat tend to drive the lapae rate on both sides of 
the front toward neutral, but the differential thermal advectkm by the ageostrophic 
wind dominates, stabilizing the warm sector and destabilizing the cold air (Ke 3 rser 
and Anthes, 1982; herdnafter, KA82). 

Cold fronts slope toward cdd air with slopes on the order of 50-100:1 (i.e., 
they ate nearly horizontal). Warm fronts have slopes toward warm air on the order 
of 200:1. The sense of the along-front temperature gradient also distinguishes warm 
fronts from cold fronts. For a cold front parallel to the y axis, dTfdy < 0 and, from 
the thermal wind relationship, du/dz > 0 while warm fronts have just the reverse 
(Keyset and Pecnick, 1987). The sense of the along front temperature gradient has 
important implications for frontogenesis. The action of cyclonic circulation on the 
along front temperature gradient is frontolytic for warm fronts and frontogenetic 
for cold fronts (Gidel, 1978; Keyser and Pecnick, 1987). Cold and warm fronts 
also often differ in terms of static stabiUty; the dr behind the cold (warm) front is 
typically unstable (stable). Levy (1989) proposes a mechanism whereby differences 
in stability result in more rapid warm frontogenesis m the early stages of frontal 
development. In some situations, the warm sector is unstable. In these cases, there 
may be strong frontogenesis (Nuss, 1989) and an intensification of the ageostrophic 
transverse vertical circulation (Koch, 1984; Reeder, 1986). 

An important aspect of frontal structure is the ageostrophic transverse ver¬ 
tical circulation about the front. The Sawyer-Eliassen equation is a linear partial 
differential equation which expresses the stream function of the ageostrophic circu¬ 
lation in terms of the geostrophic momentum and thermal fields. This ageostrophic 
stream function is consistent with a thermally direct circulation centered on the 
front (Keyser and Shapiro, 1986). A review of Sawyer’s analysis is presented by 
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Eliassen (1960) with derivations of modifications to the Sawyer-Eliassen equati<» 
for different forcings to the transverse circulation. The sense and positicm of the 
vertical circulation is sudi that the front is looked in a zone of convergence at the 
surface which helps to intensify it (Williams, 1972). Hoskins (1982) and Sh^>iro 
(1987) postulate a feedback between front<^enetic forcing and the vertical drcula- 
tion. Analyses by Sanders (1955), Ogura and Portis (1982), and Bond and Fleagle 
(1985) show that the ascent ahead of the front b concentrated into an intense ver¬ 
tical jet while the descent to the rear b more diffuse. Thb feature b apparent also 
in simulations of fronts by KA82 and Orlanski et al. (1985). 

D. FRONTOGENESIS 

1. Quasi-Geostrophic IVontogenesis 

The derivation of the quasi-geostrophic equations proceeds from the as¬ 
sumption that the Rx>ssby number (V/fL) is small. Thb b a fairly restrictive as¬ 
sumption and the equations cannot be expected to describe the evolution of small- 
scale features such as fronts. As dbcussed by Williams (1967), the conservation of 
quasi-geostrophic potential vorticity by the quasi-geostrophic system implies that a 
front at which potential vorticity is large cannot evolve from a state in which the 
potential vorticity is everywhere small. Nevertheless, Williams and Plotkin (1968), 
Williams (1968), and Stone (1966) used the quasi-geostrophic equations to inves¬ 
tigate frontogenesis forced by a horizontal deformation field. Stone treated only 
certain limiting cases. Wilhams and Plotkin used a time dependent, Boussinesq, 
adiabatic, and invisdd model with rigid surfaces at the top and bottom. They 
included an initial temperature field which contained a temperature gradient that 
was invariant with z and was confined to a zone of finite width. Thb work was 
extended by Williams (1968) to include vanishing temperature gradients at infinity 
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•ad man general boundary conditkms. In tbeei* studiee, discontinuities formed <mly 
as f oo and only at the boundaries. Away from the boundaries, the temperature 
gradient increased toward steady-state values and the scale was the Rossby radius 
of deformation. The fronts did not slope, had infinite horizontal wind speeds, and 
zero relative vorticity at the front (with la^ vortidty to either side). The lack 
of a slope in the vertical lead to regions of negative static stability. There was a 
thermally direct vertical circulation. 

2. Semi-Geostrophic iVontogenesis 

There are many situations in which the scaling arguments used in the 
development of quasi-geostrophic theory are not appropriate. In some of these situa¬ 
tions, the Rossby number is 0(1) or larger. In others, the scaling is anisotropic. One 
situation in which both of these complications are present is atmospheric frontogen- 
esis. By introducing appropriate length scales in the along-front (L) and cross-front 
(1) directions, and nondimensionalizing the equations of motion, it can be shown 
(Pedlosky, 1979; Gill, 1982) that the acceleration in the cross front direction is 
0(1/L) 1 and therefore it can be neglected. Thus, there b geostrophic balance in 

the cross-front direction. This approximation, which is referred to by Eliasen (1948) 
as the geostrophic momentum approximation, greatly simplifies the equations and 
provides for the advection of the geostrophic momentum by the total wind. 

It is convenient also to introduce the “geostrophic coordinates” (X, Z, T) == 
(z -I- v/f, z,t). For a fluid parcel in purely geostrophic motion, x = X (HB72). The 
Jacobian of this transformation is the absolute vorticity and the inverse transforma¬ 
tion- is 

z = X-vff 
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so that the absolute vorticity in (z, z,t) space is pven by 

^~ax 

Thus, when the relative vorticity is equal to / in geostrophic coordinate space, the 
absolute vorticity is infinite and the transformation is discontinuous. In this manner, 
formation of a discontinuity will occur with a relatively modest increase in vorticity 
whidi can be obtained in either the deform^ion case or the growing Eady wave 
case. This is in contrast to the quasi-geostrophic system, which will not produce a 
discontinuity at all in the Eady wave case and requires infinite time for formation 
of a discontinuity in the deformation case. 

The primitive equations incorporating the geostrophic momentum ap¬ 
proximation and transformed to geostrophic coordinate space are called the semi- 
geostrophic equations. The semi-geostrophic system is iq>plicable to a much wider 
class of phenomena than the quasi-geostrophic equations because of the less restric¬ 
tive character of the derivation (Hoskins, 1975). Since the ageostrophic velocity is 
implicit in the transformation of coordinates, the behavior of the system can be 
more easily analyzed. FVirther discussion and extensions of the semi-geostrophic 
system are given by Blumen (1981) and Schubert (1985). 

Several authors have discussed semi-geostrophic frontogenesis. Common 
to most of these studies is the use of the Boussinesq approximation, a domain 
bounded top and bottom by rigid surfaces at which the vertical velocity vanishes, 
and a constant Coriolis parameter. The discontinuity in temperature always forms 
first at the boundary. This does not produce realistic upper tropospheric fronts, but 
surface fronts are fairly faithfully simulated. It can be shown (Hoskins, 1971) that, 
when the potential vorticity is constant, analytic solutions always exist provided 
that the initial tonperature field is linear in z and an arbitrary function of z. 
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Hoskins (1971) presoited scou-geostrophic sc^utions for frcmts formed 
by confluence in a deformation vrind field. Three different initial temperature fields 
were used. The first case involved warm temperature at z » oo and cold at z 3= ->oo. 
Results showed very realistic fronts forming on the warm side at both the upper and 
lower boundaries with convergence at the surface and a thermally direct circulation. 
The second case was identical to the first except that warm air was placed at the 
origin. In this case, a situation analogous to an occlusion was obtained. The third 
case was identical to the first except that cold air was placed at the origin. A front 
formed at the upper boundary only in this case. Numerical solutions were sought 
for several modified forms of the model: 

1) if the Boussinesq approximation was not made, frontogenesis proceeded more 
slowly at the surface than at the “lid”; 

2) for slowly varying potential vorticity, frt>ntogenesi8 proceeded slightly more 
slowly at the surface but very slowly at the lid; 

3) a crude parameterization of latent heat rdease indicated that precipitation in 
the rising branch of the thermally ^rect cell ahead of the front would enhance 
surface frontogenesis; 

4) incorporating Ekman layer suction to simulate surface stress resulted in a 
weaker surface front (enhanced convergence at the surface was frontogenetic, 
but compensating divergence aloft was frontolytic). 

Upper tropospheric frontogenesis was also discussed. 

The mathematical development of the techniques used by Hoskins (1971) 
was given by HB72. They also treated the deformation wind field case and numer¬ 
ically investigated nonBoussinesq effects, the effect of variable potential vorticity. 
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latent heating, boundary layer suction, and upper tropospheric fronts. In addition, 
HB72 obtuned analytic solutions to the problem of frontogenesis forced by a grow¬ 
ing nonlinear baroclinic wave (the most unstable Eady mode). This problem was 
first investigated nummcally by Williams (1967). In this study, the domain was 
* infinite in the y direction with a uniform temperature gradient, providing an infinite 
reservoir of available potential energy. Initially, the Eady wave growth rate was 
large as it drew on the infinite source of potential energy. As the wave grew, hori¬ 
zontal shears developed much as suggested by quasi-geostrophic theory. When the 
perturbation attained finite amplitude, however, the ageostrophic velocity caused 
a distortion and further contraction of the temperature field, resulting in a realis¬ 
tic front. The analytic solution obtuned by HB72 compared quite favorably vnth 
Williams’ solution, indicating that the assumption of cross-front balance is a good 
one. 

The HB72 baroclinic wave model has been compared to observed fronts 
by Blumen (1980) and Ogura and Portis (1982). Blumen compared the model to 
Sanders (1955) analysis and Ogura and Portb compared it to their analysis of a 
front from the severe environmental storms and mesoscale experiment (SESAME). 
Although the model was capable of reproducing many features of observed fronts, 
there were some unrealistic features: 

1) the model surface wind speed (50 m s“*) was too large; 

2) the model vorticity was too small; 

3) Sanders (1955) showed that the observed front was maintained as a balance 
between frontogenesis and frontolysis while the model was predominantly fron- 
togenetic; 
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4) the model vertical circulation was too weak; the ascent was less than 20% 
of the value observed by Sanders. 

The major failing was that the model did not produce an intense upward vertical jet 
ahead of the front. In an effort to remedy this situation, Blumen (1980) introduced 
an Ekman layer into the model. This enhanced the vertical velocity, but the small 
vertical gradients in the model prevented a significant contribution to frontogene- 
sis. Blumen and Wu (1983) extended this work by including a more sophisticated 
boundary layer with a specified vertical velocity at the surface. They noted that the 
vertical velocity at the top of the layer has both an inviscid component associated 
with ageostrophic divergence at that level and a component due to Ekman la}^r 
dynamics. The resultant vertical velocity distribution was more realistic, although 
the magnitude of the upward motion was only slightly larger (6 x 10~* m s“^ vs 
4 X 10“^ m s"*). 

Reviews of HB72 have been given by Pedlosky (1979), Holton (1979), 
and Hoskins (1982). The review by Hoskins is quite extensive and includes both 
two and three dimensional semigeostrophic frontogenesis and a brief discussion of 
quasi-geostrophic frontogenesis, which points out the important distinctions between 
semi-geostrophy and quasi-geostrophy. For the case of a growing Eady wave in the 
three dimensional semigeostrophic equations, it is necessary to slightly modify the 
Eady wave and introduce buoyancy gradients in the initial conditions to obtain 
realistic three dimensional fronts. 

Numerical solutions to the three-dimensional semi-geostrophic equations 
for several different growing nonlinear baroclinic waves (including the Eady case) 
were obtained by Hoskins and West (1979). With a strong westerly jet, realistic 
cold and warm frontal zones formed in the growing waves at the leading edge of the 
cold air. 
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S. Quasi-geostrophic Versus Semi-geostrophic fVontogenesis 

It is instructive to compare the approaches taken in the previous two sec¬ 
tions. Comparison of quasi-geostrophic and semi-geostrophic frontogenesis foUowing 
Orlanski et al. (1985) serves to highlight three important dynamical processes which 
are strongly modulated by boundary layer physics. These processes are 1) a feed¬ 
back between increasing frontal gradients and increasing ageostrophic convergence, 
2) increasing convergence leads to growth in vorticity, and 3) a braking mechanism 
is postulated based on an imbalance between the mass and wind fields in frontoge¬ 
nesis. Expressions for the substantial time derivatives of the cross front potential 
temperature gradient and vertical wind shear can be written 



where an asterisk (’^) indicates semi-ge(»trophic terms added to the basic quasi- 
geostrophic system and two asterisks (4>4>) indicate terms of the full nongeostrophic 
system. 

The fundamental difference between quasi-geostrophic and semi-geostrophic 
frontogenesis is that the sum of the geostrophic and the ageostrophic wind is used 
in the horizontal and vertical advection of temperature and geostrophic momentum. 
This is apparent from the third and fourth terms on the RHS of equations (2.1) and 
(2.2). This implies that, as the frontal gradients increase, ageostrophic convergence 
also increases, which brings about stiU stronger frontal gradients. Evidently, this 
process will bring about a discontinuity in a rdatively short time. If only the first 
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two tenns on the RHS are retained (the ({uan-geoatrophic case), this ptooets cannot 
occur. The vwticity equation can be written as 

i(C) = -/C-CO‘-(V.-VCr (2.3) 

where D is the divergence. The second term on the RHS implies that cyclonic vor- 
tidty will increase in areas of convergence and decrease in areas of divergence. Thus, 
this is another mechanism which will enhance ^ontogenesis in the semi-geostrophic 
system relative to the quasi-geostrophic qrstem. 

A shortcoming of the semi-geostrophic system is that the continual action 
of the vortex stretching term (the second term on the RHS) leads to unbounded 
growth in vorticity. Use of the full nongeostrophic equations provides a possible 
braking mechanism. The <hvergence equation is 

(2.4) 

The imbalance between / and p (which is essentially an imbalance between the mass 
and the wind fields due to ageostrophic motion *) may resiUt in a positive contribu* 
tion to divergence or a decrease in convergence. This may reduce firontogenesis and 
(from terms 1 and 2 on the RHS of 2.3) slow the production of cyclonic vorticity. 
In a high resolution nonhydrostatic simulation of frontogenesis. Gall et al. (1987) 
found some indication of this braking mechanism (manifested by a separation of 
the re^on of maximum convergence from the region of maximum vorticity). How¬ 
ever, the cyclonic vorticity continued to increase. Using a simple nonhydrostatic 
model. Levy and Bretherton (1987) critically examined the braking mechanism and 
found that, in several different plausible situations, the medianism was frontogenetic 
rather than frontolytic. They postulated that the phase shift between the region 
*Ia quast-geortrophic and aeini-feoatKq>hic frcwtagencais, a balance exiatt between /C and V^p. 
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of maximum convergence and the region of maximum vortidty might be due to 
gravity waves excited by frontogenesis. In a comment on this paper, Gamer (1989) 
criticized the author’s analysis as being irrelevant and maintained that the gravity 
waves are too weak to bring about the observed phase shift. Levy and Bretherton’s 
reply is not compelling. In a two dimensional investigation of warm and cold fronts, 
Keyser and Pecnick (1987) found no evidence of a dynamical brake on frontogenesis. 
Thus, further study will be required to verify this braking mechanism. 

4. FVontogenesis With the Primitive Equations 

Several authors have examined frontogenesis using the primitive equa¬ 
tions. These studies can be viewed in a hierarchy of increasing sophistication. 
Inviscid, adiabatic, Boussinesq, hydrostatic models were used by Williams (1967) 
(discussed above) and Williams (1972). Williams (1972) used a horizontal defor¬ 
mation field in a comparison of quasi-geostrophic and nongeostrophic frontogenesis, 
essentially complementing earlier work (Williams, 1968). Results indicated that the 
nongeostrophic nonlinear approach provided more realistic fronts with vertical tilt 
and the tendency to form a discontinuity in finite time. Parameterizations of hori¬ 
zontal and vertical diffusion of heat and momentum were added by Williams (1974), 
relaxing the adiabatic and inviscid condition. These parameterizations included 
only the neutral case, involving constant eddy coefficients and a simple Ekman 
layer. Inclusion of surface stress and turbulent diffusion of heat in this study re¬ 
sulted in the production of quasi-steady state fronts after 1-2 days of integration. 
Vertical diffusion of heat was found to play a dominant role in determining the 
steady state structure. NonBoussinesq (fully compressible) effects were introduced 
by K A82 along with more elaborate boundary layer physics. A two dimensional ver¬ 
sion of the Pennsylvania State University/National Center for Atmospheric Research 
(PSU/NCAR) model was used to investigate the effects of differing boundary layer 
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parameterisAtions on frontogenesu. A bulk boundary layer and a hi^ reMluti<»i K- 
theory boundary layer were used. The high resolution boundary layer provided the 
sharpest front. This model was also used by Kuo et al. (1991a) with the high reso¬ 
lution boundary layer. These studies will be discussed in more detail below. Finally, 
the nonhydrostatic primitive equations were used by Gall et al. (1987) in extremely 
high horizontal and vertical resolution inviscid integrations designed to determine 
the minimum scale of fronts. Results indicated that continuously decreasing the 
vertical mesh size resulted in smaller and smaller frontal scales (results were rela¬ 
tively insensitive to horizontal mesh size given the shallow slope of the front). Thus, 
no minimum scale could be determined. The introduction of horizontal and vertical 
diffusion would modify this conclusion. Inviscid nonhydrostatic primitive equation 
simulations were also performed by Polovarapu and Peltier (1990) with much lower 
horizontal and vertical resolution. They investigated differences between simulations 
of baroclinic wave life cycles using the /-plane and ^-plane approximations. 

E. IMPORTANCE OF BOUNDARY LAYER PROCESSES 

In general, one can envision that surface stress, static stability, and turbulent 
diffusion must be important in the formation of zones of large horizontal gradients 
in the atmosphere. In particular, the formation of discontinuities is not observed 
in models incorporating turbulent diffusion since turbulent mixing will reduce the 
magnitude of gradients in temperature and wind speed. Hoskins (1971) and HB72 
note that the existence of a Richardson number near the critical value in the vicin¬ 
ity of modeled fronts implies that turbulent mixing will be important. Steady-state 
fronts were produced by Williams (1974) using simple parameterizations of horizon¬ 
tal and vertical turbulent diffusion of heat and momentum. Using the PSU/NCAR 
model in a study of explosive marine cyclogenesis, Kuo et al. (1991b) noted that 
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the structtire of modeled fronts was a great deal more realistic when boundary laya 
physics was included. The physical parametoizations involved the Blackadar (1979) 
K-theory scheme and surface sensible and latent heat fluxes based on similarity the¬ 
ory. In reality, fronts exist in quasi-steady state for portions of their lifetime due 
to a balance between the frontogenetic effects of large-scale convergence and the 
frontolytic effects of the boundary layer. 

One of the most important features of frontal structure is the thermally direct 
ageostrophic vertical circulation in the cross-front plane.* Convergence associated 
with this circulation intensifies cross-front gradients and it increases in magnitude 
as they do. Several observational and numerical studies have shown this circulation 
to be forced by the boundary layer. Carbone (1982) demonstrated (using Doppler 
radar and surface-based observing systems) that the strong vertical updraft in fronts 
observed in central California resulted from boundary layer forcing. In his compar¬ 
ison of Sanders (1955) analysis to the HB72 model, Blumen (1980) introduced an 
Ekman layer into the model in an effort to improve the vertical velocity field. Blu¬ 
men and Wu (1983) further modified the model by specifying the vertical motion 
at the top of the boundary layer at the initial time in an effort to bring the verti¬ 
cal velocity closer to the observations. They concluded that the low-level vertical 
motion field was very sensitive to the boundary layer parameterization used. Bond 
and Fleagle (1985) found boundary layer convergence to be a major component in 
forcing vertical motion in fronts and storms over the NE Pacific Ocean. Moreover, 
Fleagle et al. (1988) found that 80% of the upward vertical motion in a cold front 
over the NE Pacific ocean was attributable to surface forcing. 

In their investigation of the effects of impact of boundary layer physical pro¬ 
cesses on frontogenesis, KA82 utilized a two dimensional version of the PSU/NCAR 
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modd with » frowing nonlinear barodink wa.ve in the initial coa<htk»a and per¬ 
formed four experimoits: 

1} adiabatic and inviacid with moderate vertical resolution^ 

2) bulk-drag boundary layer with moderate vertical resolution; 

3) adiabatic and inviscid with high vertical resolution (approximately 5 levels in 
the lowest km); 

4) high resolution K-theory boundary layer and thermally insulated lower bound¬ 
ary with high vertical resolution. 

They did not consider the effects of moisture, surface heat flux, ot radiaticm. Ra¬ 
diation was considered to be of no importance in frontogenesis by HB72. Hoskins 
(1982) stated that, while moisture may be important in the structure of steady-state 
fronts, it b not necessary to include it in studying frontogenesb. 

The results of KA82 indicated that adding boundary layer physics as in 4) 
above to the adiabatic and inviscid simulation resulted in almost immediate produc¬ 
tion of large vertical shears in u (cross-front) and v (along-front) whidi caused vig¬ 
orous mixing. FVictional depletion of v momentum resulted in a large ageostrophic 
V component and large cross isobaric acceleration in u, which invalidated the as¬ 
sumption of cross-firont balance. A quasi-steady state devebped after about 8 h. 
Vertical fluxes of sensible heat at levels above the surface tended to drive the l^we 
rates on both sides of the front toward neutral but vertical differential thermal ad- 
vection dominated. The frictkmally driven ageostrophic boundary layer inflow was 
responsible for contraction of the temperatiire gradient and thus counteracted the 
dissipation of potential vortidty due to friction. The frictionally driven ageostrophic 
boundary layer inflow was also responsible for the creation of a vertical jet ahead 



of the surface front after about 8 hours. Thus, the high resolution boundary layer 
formulation provided a mudi more realistic vertical vdodty fidd than the bulk 
boimdary layer. 

Although the study of KA82 succesfully reproduced many features of observed 
• frontogenesis, it suffers fron several drawbadcs. Most of these ston from the use 
of the insulated surface condition. For the bulk formulation, the maximum heat 
flux is required to be at the surface. It is not clear how this can occur with an 
insulated surface; the bulk formulation applks only in convectively driven, well- 
mixed conditions. For the high resolution boundary layer, the heat flux at the 
surface was assumed to be zero and only the nocturnal case was used due to the 
insulated surface. This is clearly inappropriate in the well mixed area bdund the 
front. For these reasons, the interpretations given by the authors, specifically those 
regarding the evolution of the thermal field, may not be justified. In particular, the 
dominance of vertical differential thermal advection over vertical fluxes of sensible 
heat in controlling static stability requires further investigation. The authors state 
in their conclusion that three dimensional simulations will eventually be required to 
answer some of the questions rdsed in thb work. The addition of a third dimension 
will allow more realistic representation of horizontal divergence which may alter the 
vertical velocity field. Finally, the authors do not discuss in any detail the differences 
in model frontogenesis caused by the differing boundary layer parameterizations. 

Although there have subsequently been a number of similar studies, none has 
been as broad in scope at that of KA82. Several other authors have investigated 
the effect of surface sensible heat fluxes on frontogenesis. Moore (1991) included 
a surface heat flux parameterization into the two dimensional HB72 framework. 
Results indicated that the addition of surface heat flux reduced the stability in 
the frontal zone and, for certain surface temperature distributions, increased frontal 


19 






• III . 

4 


buoclinicity. This study omits prooessa acoouated for in the foil primitive equatioo 
framework and three dimensional effects. The impact surface drag on frontogen- 
esis was investigated by Hines and Mechoso (1993). They used a constant surface 
drag coefficient in an otherwise invisdd simulation of frontogenesis forced by a grow¬ 
ing baroclinic wave. Thrir results showed that a feedback between deformation and 
warm frontogenesis was sensitive to surface drag. With larger drag coefficients, the 
feedbadc was inhibited and the warm front weakened. In an idealized simulation of 
a warm front near the Gulf Stream, Nuss (1989) found that the effect of stability on 
surface stress was critical. In fact, he found that a stratification gradient across the 
front was necessary to enhance frictional convergence and frontogenesis. Stability 
was also found to play a significant role in the investigation of the pre-ERICA storm 
by Neiman et al. (1990). In this case, the front moved past the Gulf Stream over 
the cold ('^ 0 C) waters of the North Atlantic. Stable stratification in the warm 
sector resulted in decoupling of the momentum from the surface, giving rise to a 
strong low level jet (~ 40 m s~^ at 'S' 500 m). Thermal advection associated with 
the LU strengthened the fronts above the surface. In light of the results of Nuss 
and Neiman et al., the use of a constant drag coefficient by Hines and Mechoso 
seems to be an unwarranted simplification. In the investigation by Nuss, the sea 
surface temperature distribution was specified so as to represent the large thermal 
gradients typical of the Gulf Stream. The generality of these results is thus unclear. 

Thus, a numbo* of previous studies have examined various aspects of the im¬ 
pact of boundary layer processes on frontogenesis and the thermally direct ageostrophic 
vertical circulation. None of these studies, however, have addressed the issue of the 
specific processes which must be correctly parameterized in order to c^ture the 
essential details of maritime atmospheric frontogenesis in three dimensions. Only 
the two-dimensional study by KA82 addrened the impact of differing parameteri- 
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zations of the boundary layer on frontc^enem. This study was flawed by the use 
of an insulated lower boundary. Other investigations of the impact of surface fluxes 
have been inconclusive in that the models used were highly idealized and therefore 
unable to assess the nonlinear interactions of surface heat flux with other physical 
processes. These considerations lead to the basic hypothesu to be examined in the 
present study. 
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III. DESCRIPTION OF PROPOSED WORK 

The central hypothesis of this study is that increasing the sophistication of 
the boundary layer parameterization will increase the ability of the model to repro¬ 
duce the three dimensional features of observed frontogenesis. In particular, it is 
hypothesized that more realistic representation of the boundary layer will provide 
an improved simulation of the ageostrophic vertical circulation and frontogenetic 
processes associated with it. 

In an attempt to verify this hypothesis, three different versions of the Navy Op¬ 
erational Regional Atmospheric Prediction System (NORAPS) will be used. These 
versions will consist of an adiabatic and invisdd run and two turbulence parame- 
terizations of differing levels of sophistication. The model is a hydrostatic primitive 
equation model with a scheme C grid in the horizontal and a a(s p/p«) coordinate 
system in the vertical. Split explicit time differencing (in which the fast gravity 
modes are treated separately) is also used. As in KA82, radiation and clouds will 
not be included. Radiation is of little importance and, while moisture may play an 
important role, the intent of the present study is investigation of boundary layer 
turbulent processes occuring near the surface; inclusion of clouds could complicate 
interpretation of the results. Omission of clouds and radiation will also significantly 
decrease computational expense. The initial conditions will consist of a highly baro- 
clinically unstable state. The initial wind, temperature, and geopotential fields will 
be in hydrostatic and thermal wind balance. Prior to the start of the computational 
cycle, a small amplitude perturbation will be superimposed on the v wind compo¬ 
nent and geopotential fields. The lower boundary conditions will be specified so as 
to simulate an oceanic surface. 
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vcfsions of the model are: 

1) Adiabatic and inviscid 

This version will serve as an experimental control for comparison to versions with 
physical processes active. 

2) K-theory 

In this version, the turbulent fluxes of heat and momentum ate parameterized by 
assuming that the fluxes are proportional to the gradients of the mean quantities 
(with K being the coefficient of proportionality). The K>theory version involves a 
mixing length based on Blackadar (1962) whidi asymptotically approadies a con¬ 
stant value and the eddy coefficients are based on the functional dependoice on the 
Richardson number devised by Louis (1979). This version is similar to that used 
by KA82 in that the value of K is a function of the local Ridiardson number. This 
experiment will reveal the impact of surface fluxes on a simulation similar to that 
performed by KA82. The ctixrent investigation represents an extension of the work 
of KA82 also in that three-dimensional aspects of the problem are addressed. 

3) Second-order closure boundary layer physics 

The second-order closure scheme makes explidt use of the ensemble average vari¬ 
ances and covariances of the turbulent heat and momentum fields (turbulent fluxes). 
This version uses the Mellor and Yamada (1974) level 3 model, in which prognostic 
equations are used for turbulent kinetic energy and the mean quantities as well as for 
the variance and covariance of the turbulent temperature and moisture fluctuations. 
Additional diagnostic equations are induded for the turbulent mom^tum fluxes. 
This will provide for a furly realistic treatment of the turbulence, particularly when 
the turbulence level is changing rapidly. Rapid changes in turbulent kinetic energy 
are to be expected in frontogenesis; as the strength of the front increases, the wind 
shear will also increase, as will the buoyant production of TKE. Thus, both the 
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mechanical and buoyant production tmns of the TKE equation will be large, lead¬ 
ing to a large local tendency in TKE. The conservative variables are liquid water 
potential temperature ^ (in this case, liquid water potential tempaature reduces to 
ordinary potential temperature), total water substance, u momentum, and v mo¬ 
mentum . The prognostic equations are solved semi-iiiq>licitly using the Thomas 
algorithm (Richtm^er and Morton, 1967). This experiment will serve to demon¬ 
strate that additional insights can be obtained using a turbulence parameterization 
more sophisticated than that used by KA82. 

The sensitivity to rapid temporal variations in turbulence level will be investi¬ 
gated by using a level 2 version of the model. In a Yamada and Mellor (1982) level 
2 scheme, prognostic equations are solved for the mean quantities, but diagnostic 
equations are used for turbulent kinetic energy and all other turbulence quantities 
such as the variances of temperature and moisture and the temperature/moisture 
covariance. 

A detailed analysis of frontogenetical forcing will be performed for each of 
the different parameterizations. It will be possible to determine what processes are 
most important in frontogenesis (and frontolysis) for each of the parameterizations 
at various stages in the development of the front. 

In order to determine which version gives the most faithful representation of at¬ 
mospheric frontogenesis, the results can be qualitatively compared to observational 
studies. The analysis of frontogenetical forcing will facilitate semi-quantitative com¬ 
parisons with atmospheric fronts; similar analyses have been performed for observed 
fronts by Sanders (1955), Ogura and Portis (1982), Carbone (1982), and Shapiro 
(1984). 

^$ 1 ^ = 0 — {O/T^L/CfQi, where Qi is liquid water specific humidity. 
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Sensitivity studies will also be conducted in <»der to assess the unpwtanoe 
of surface heat and momentum fluxes. While KA82 did not include the effect of 
surface flusms, it has been suggested recently (Dorian et al., 1989; Moore, 1991) that 
horizontally differential surface sensible heating may be important in firontogenesis. 

^ None of these studies, however, provide a complete understanding of the impact of 
surface fluxes in the maritime case. The importance of surface drag is addressed by 
Nuss (1989) and Hines and Mechoso (1993). In the present study, a more realistic 
distribution of surface drag will be used than in Hines and Mechoso and the sea 
surface temperature field will be more typical of maritime firontogenesis in the open 
ocean than that used by Nuss. 

The results of this study should provide important insights into the three di¬ 
mensional nature of the boundary layer processes important in frontogenesis. The 
evolution of the vertical updraft forced by the inflow with model-generated hor¬ 
izontal divergence can also be examined. The frontolytical action of diffusion in 
combination with frontogenetic aspects of the vertical updraft and vertical gradi¬ 
ents above the boundary layer can be studied. The sensitivity study of the surface 
momentum and sensible heat fluxes should shed some light on these features as well. 
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IV. NUMERICAL EXPERIMENTS 


A. REGIONAL FORECAST MODEL 

The forecast model is a hydrostatic primitive equation model with <r(= p/x; x = 
surface pressure) as the vertical coordinate. The grid is staggered following sdieme 
C. This C staggering involves placing the U component of the wind one half grid 
interval in the x-direction away from mass points (the points at which teiiq>erature, 
pressure, moisture, and vertical velocity are defined) and the V component of the 
wind one-half grid interval in the y direction away from mass points. In the vertical, 
the vertical velocity dajdt is defined on layer interfaces and all other quantities are 
defined at the center of each layer (Arakawa and Lamb, 1977). Time integration is 
performed using the split-explicit method. In this technique, prognostic equations 
are solved with a centered time difference with forcing due to the meteorological 
modes and slowest gravity modes. These solutions are then adjusted for the linear 
terms governing the fast gravity modes (Madala, 1982). Robert time filtering is ap¬ 
plied at the end of each time step to prevent solution splitting. Robert time filtering 
is described by Haltiner and Williams (1980). Fourth-order advection and diffusion 
are used for all of the prognostic variables (second-order diffusion is used on the 
first interior row). Diffusion computations on sigma surfaces will lead to spurious 
sources and sinks of energy when applied to mass and moisture fields. To overcome 
this problem, the fourth-order diffusion operator is applied to deviations from the 
standard atmosphere for temperature and moisture. Further details concerning the 
model can be found in Hodur (1987). 
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The primitive equetions can be e xpre w e d in or ooordinatei u 
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dt 

d{xuv,) d{xw>) d{xua) 

dx dy da 
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(4.1) 
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(4.3) 

dixq) d{xuq) dixvq) dixqa) 

dt ^ dx dy da ^ 



+Tjrvj(vj(, -,,)) 

(4.4) 



(4.5) 



(4.6) 



(4.7) 


Equations (4.1) and (4.2) are the horizontal equations of motion, equation (4.3) is 
a statement of the first law of thermodynamics, equation (4.4) is the conservation 
equation for water substance, equation (4.5) b the hydrostatic equation, equation 
(4.6) b the continuity equation, and (4.7) is the pressure tendency equation. 

Characteristics of the model and initialization procedure are summarized in 
Table 4.1 and the distribution of computational levels in the vertical is given in Table 
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4.2. With the exception of the a vertical coordinate, the first six iteoas in Table 
4.1 were chosen based on a balance between a desire to resolve the phenomena 
under study and acceptable running time and memory allocation. The jet speed 
and perturbation amplitude were chosen so as to pve a convmient baroclinic wava 
growth rate. The horizontal diffusion co^cient is quite conservative and results 
in very little diffusion. The distribution of votical levels in Table 4.2 was chosen 
so as to give high vertical resolution near the surface and near the tropopause; 
high vertical resolution in the boundary layer is crucial in this study and failure to 
adequately resolve the tropopause can lead to numerical difficulties. 

B. INITIAL AND BOUNDARY CONDITIONS 

The initial conditions for the numerical experiments to be performed consist of 
a strongly baroclinically unstable regime. This regime involves a 40 m s~^ westerly 
jet at 200 mb at the center of the domain, decreasing to zero at the north and 
south boundaries with a sine^ profile. In the vertical, the u component of the wind 
decreases (linearly in z) to zero at the surface. The jet, which is invariant in the 
zonal direction, is in geostrophic and thermal wind balance. This implies a strong 
meridional temperature gradient at all levels with a concentrated baroclinic zone 
near the center of the grid. The v component consists of an altitude independent 
3000 km sinusoidal perturbation in both x and y with a maximum amplitude of 1 
m s~^ at the center of the domain. The wavelength of the perturbation was chosen 
to be near the wavelength of maximum instability from linear theory (3246.98 km). 
The vertical temperature profile was chosen so as to yield a value of the stability 
parameter which enhances the growth rate of a 3000 km baroclinically unstable 
wave. The temperature and pressure gradients and both velocity components are 
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800 «t the noth and south boundaries and all fields are periodic in the east>west 
direction. 

The initialization procedure involves first qtedfying the jet profile at eadi levri. 
Then, using the velocity field at the lowest levd, a modified form of the balance 
equation is applied at the surface to obtain surface pressure. Using the surface 
pressure and wind field, the balance equation is iq;}plied at each level to obtain 
the geopotential field. Finally, an altitude independent, sinusoidal perturbation is 
added to the geopotential field and the geostrophic perturbation v component is 
determined. 

As in Eady’s (1949) investigation of baroclinic instability, the model equations 
are formulated on an /-plane = 0) and the vertical motion on sigma surfaces 
(da/di) vanishes at both the upper and lower boundaries. This implies that u; 0 
at the surface and cj = 0 at the upper boundary. 


C. BOUNDARY LAYER PARAMETERIZATIONS 


1. Introduction 

The terms in equations (4.1)-(4.6) involving F, Qy and Qm represent vis¬ 
cous and diabatic effects of the boundary layer, radiation, and stable and convective 
precipitation. The purpose of the boundary layer parameterization is to quantify 
these terms by evaluating the divergence of the turbulent fluxes of heat, moisture, 
and momentum from the surface. The rate of change of any prognostic variable i 
due to the boundary layer physics is 


dt dz 


(4.8) 


where {4>) indicates an ensemble mean and denotes a turbulent deviation from 
the mean. Viewed in this context, the boundary layer parameterizations used in the 
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presmt study represent differing levds in tl^ sophistication of the tediniques used 
to evaluate equation (4.8). 

2. Adiabatic and Nearly Inviscid Case 

The adiabatic and nearly inviscid case involves no boundary layer par 
rameterization (it is nearly inviscid in the sense that the fourth-order horizcmtal 
diffusion operator in equations (4.1)-(4.4) is included for oonq>utational stability; 
this pseudo-viscous effect is always present). The lower boundary is insulated and 
free-slip. No exchanges of heat, moisture, or momentum take place at the surface. 
This case serves as a control experiment for this investigation. 

A convective adjustment scheme is used in order to ensure computational 
stability. It can be shown that the existence of a super adiabatic lapse rate in 
any column will lead to an exponential growth in the vertical velocity. Convective 
adjustment, applied at the end of each timestep, removes super adiabatic lapse rates 
in such a way as to conserve the total potential energy. The problem of changing 
the temperatures to produce a neutral lapse rate subject to the constraint of total 
potential energy conservation reduces to a system of linear equations which is solved 
using a standard matrix inversion technique (Carnahan et al., 1969). 

3. K-Theory Boundary Layer Parameterization 

In the K-theory boundary layer parameterization, the boundary layer is 
explicitly resolved and the surface fluxes are based on Monin-Obukhov similarity 
theory. No a priori assumption about the flux profile (above the surface layer) is 
required. 

The surface layer similarity approach used to obtain the surface fluxes 
was developed by Louis (1979). This method uses analytic functions of Ria of the 
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ul^a*U*F{RiB) 


(4.10) 


11.9. * a/RUA$F{RiB) 


(4.11) 


where 


and a is analogous to Cos *nd 6, 6, c*, and R are constants. 

The surface roughness, zot u determined from Chamock’s formula for 
flow over water: 

*0 = culfg (4.12) 

where c = 0.05. 

K-theory owes its origin to a physical analogy between turbulent diffusion 
in the boundary layer and molecular diffusion. In both processes, the flux of a 
quantity ( is proportional to the gradient of that quantity; i.e., 

(4.13) 

Where K is the eddy diffusivity. The basis of the assumption that the turbulent 
flux is proportional to the gradient lies in the mixing length hypothesis. 

In the present study, the form for K is that proposed by Louis (1979): 

AV 


(4.14) 












where the mixiiig length t is givra by BUdceder (1962) 
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l + fc*/A 


(4.15) 


and 


Bxb^ 


g Ax A0 
0{AV)* 


(4.16) 


The fiuncti<m F in equatkm (4.14) is the same as in the surface layer (equations 
(4.10) and (4.11)) except that c is pvai by 


c = ((I + Ai/2)''’ -1)*^ 

Note that the mixing length i in equation (4.15) is proportional to 
height near the surface but asymptotically approadies a reference value A at higher 
elevations: 

kx \ 

1 + ibs/A " A/*z+l’ 

for kz X,t kx 

for kx > Xyi A 

This behavior is consistent with physical reasoning; near the surface, the size of 
eddies is constrained by proximity to the surface while aloft, the mixing length does 
not increase without bound with elevation. In Louis (1979), A was fpven a constant 
value of 100 m. Experience with this version of the parameterization, however, has 
shown that values of A = 150 m for Km and A = 350 m for Km yield satisfactory 
model performance for a wide variety of meteorological situations. Use of these 
values in simulating frontogenesis is not inappropriate. 

Having specified the form of /if, equation (4.8) can be solved. 

No convective adjustment scheme is necessary since the parameterized 
turbulent mixing will quickly respond to the presence of super-adiabatic lapse rates 
and re-distribute heat accordingly. 
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4. Second-Order Cloeure Boundary Layer Parameterisation 

In this vertiony no a priori assumption about the rdationship between the 
fluxes and vertical gradients is required; the variances and covariances of momentum, 
heat, and moisture are used directly in detmnining the boundary layer structure and 
evolution. In order to close the system of equations, expressions must be obt^ed 
for the fluxes. These expressions are derived from the turbulent equation of motion 
and scalar conservation, which, in turn, are derived from the conventional equations 
by introducing the decomposition 

^ = C + (4.17) 

and applying Reynolds averaging. Expressions for the fluxes are then obtained 
by cross multiplying the turbulence equations by turbulent momentum and scalar 
property fluctuation. The resulting expressions are 



The difficulty with this approach is that we have now obtained expressions for the 
second-order fluxes UjUj and in terms of the higher order moments Were 
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we to continue in this vein to derive expressions for third-<xder tenns, the resulting 
expressions would contiun fourth-order tarns. This is the fundamental nature oi the 
turbulence closiure problem. Once the turbulent field is decomposed into a mean 
and a perturbation (using equation (4.17)), the full details concerning the flow’s 
instantaneous behavior can only be reconstructed with knowledge of all ensemble- 
average statistics of all turbulent moments (Burk, 1978). In second-orda closure, 
we retain only the first two turbulent moments and parameterize the third-order 
moments in terms of the second-order moments, thus obtaining closure. This is in 
contrast to K-theory in which the second-order fluxes are parameterized in terms of 
the mean quantities (K-theoiy is sometimes referred to as ‘iirst-order closure”). 

The second-order closure boundary layer parameterization utilizes a ver¬ 
tically nested grid structure in which the physics (turbulent diffusion, radiation, and 
moist thermodynamics) computations are performed on a high vertical resolution 
grid nested within the regional model grid. The dynamics (horizontal and vertical 
advection and coriolis and pressure gradient) computations are performed on the low 
vertical resolution regional model grid. In the present study, the vertical nesting is 
not used. While the code used to perform the interpolations from the fine grid to 
the regional model grid is active, there exists a one-to-one correspondence between 
the fine grid and the regional model grid. This simplification has been introduced 
in order to facilitate comparison with the K-theory simulation. 

The surface layer similarity approach used to obtain the surface fluxes 
was developed by Liu et al. (1979). This approach implicitly incorporates an inter¬ 
facial sublayer or inner surface layer in which molecular diffusion is important. It 
addresses certain fundamental differences between the way in which scalar quantities 
(e.g., heat and water vapor) and momentum are transported in this inner surface 
layer. The bulk transfer coefficients for scalar quantities are vastly different from 
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those for momentum. Thk difference it, in part, due to different rougjhnett lengtht 
computed for mcuneotum (re) and tor heat and water vi^mt (rt, zq). 

The differences between the transport of momentum and that of heat 
and water vi^r are primarily due to the action of pressure forces (normal stresses) 
on the roughness elements which transport momentum. Thus, momentum is trans¬ 
ported in the inner surface layer by both molecular diffusion and pressure forces 
while heat and water vi 4 >or are transported only by molecular diffusion. This re¬ 
sults in a decrease in the bulk transfer coeffidents for heat and water vapor with 
increasing wind speed above a threshold value. The transfer coeffident for nmmen- 
tum increases monotonically with increasing wind speed. These differences in the 
transfer coefficients are due to differences in the roughness lengths, which are ob¬ 
tained by matching profiles in the outer surface layer (where molecular diffusion is 
not important) to those from the inner surface layer. 

The proiiies for the outer surface layer are given by 


(T - Z) ^ ln(z/zT)-^ 

T, ~~ auk 

{Q - Q,) _ ho) - 00 

Qm ask 

{U - U.) ^ ln(z/zo) - th 

U, k 


(4.20) 

(4.21) 

(4.22) 


where an = KhIKm and as = KeIKm at neutral stability {Kh^ Km, and Ke 
are the turbulent diffusivities for heat, momentum, and moisture, respectively). 
Matching these expressions with corresponding expressions for the inner surface 
layer yields expressions for the surface roughness: 


zx = (4.23) 

z, = (4.24) 
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Zo = Giv/u, iu,lv<rx 

*6 * ri < iu,/v < rt (4.25) 

2 o = r 2 < iu,/v 

where ai, a^, bi, 62 , Gi, G 3 are constants^ v is the kinematic viscosity, ( is the 
diameter of roughness elements, and Rr is the roughness Rqmolds number {R, as 

ZoU,/l/). 

Mellor and Yamada (1974) devised a hierarchy of second-order doeure 
modeling techniques based on the degree of turbulent anisotropy admitted by the 
equations used. In this work, which represented a landmark in second-order closure 
theory, four levels were discussed, ranging from level one in which ail quantities 
are treated diagnostically (completely algebraic) to level four in which all quantities 
are expressed prognostically (completely differential). Subsequently, Mellor and Ya¬ 
mada (1982) identified the “level 2.5’* model in which only the mean quantities and 
the turbulent kinetic energy are treated prognostically. In terms of this turbulence 
closure hierarchy, we have used both level 2.5 and level 3 closures in the vertically 
nested model. At level 3, we solve the following prognostic equations for turbulent 
kinetic energy, temperature variance, moisture variance, and temperature-moisture 
covariance: 


d /, Sen ^-r-,8V 


+20gu/6\, - 2 e®/Ai 

(4.26) 


(4.27) 


(4.28) 

T = 1 (^ 7 ’) -^'7 -^7 

(4.29) 
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(•l^caiioa of tlie boundary layer approadmatiwi shows hcNriaontal derivatives to 
be negligible in comparison to vertical derivatives). In equatuHis (4.26)-(4.29), 
e’/2 is the turbulent kinetic energy, $i is the liquid water potential temperature 
(9iss0 — $L/TCpqi)^ and is the total water substance specific humidity. At level 
2.5, equations (4.27)-(4.29) are replaced with diagnostic expressions: 


2 ^ A 



A Kk^q^dq^ 
e dx dx 


tstt k 


(4.30) 

(4.31) 

(4.32) 


The covariances which must be evaluated to solve equations (4.26)-(4.29) are u'u/, 
t/ti/, ti/dj, and w*ql,. In equations (4.26)-(4.32), L, Ai, and Aj are scaling 
factors proportional to a master length scale f. We use the closure constants from 


Mellor and Yamada (1982): 


I * 0.2f; Ai = 16.W; As = lO.lf 


Yamada and Mellor (1975) show that the equations for the momentum 
and heat fluxes may be written in the form 




dz 


_ dV 

-x/w' = 

_ aa _ 


where 


r = 


ZAite^g 




(4.33) 

(4.34) 

(4.35) 

(4.36) 


H«e, P is the coefficient of thermal expansion (= 1/^), while A\ = 0.92 and A% = 
0.74 are closure constants from Mellor and Yamada (1982). The eddy coefficients 
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Km And Kk are given in Yanutda and Mellor (1975) as algebraic functions involving 
both predicted mean and turbulence variables. 

An equation for the total moisture flux may be written as 




(4.37) 

where 

AT. = 2^53 
e 

(4.38) 

and 

„ Zaii0g 

1 ^ — 

e 

(4.39) 

The vertical velocity variance in equation (4.38) can be evaluated as 



c \dz J \az J j e 

(4.40) 


The remaining unknown quantity in equations (4.26)-(4.29) is 
Yamada and Mellor (1979) show that 

(4*41) 


where 

/^ = 1 + .6099w — 1.609gj 

V BL 

01 = (1 + .6099,,, ~ 3.22gi)^;^— .6O991 

OL 

0^ = .609(0, + — 9 ,) (4.42) 

Equation (4.41) can be re-arranged to obtmn an equation for w’d/. How¬ 
ever, new unknowns have been introduced in equations (4.37) (0J,9i,) and (4.41) 
(u7^). Yamada and Mellor (1975) give the following expression for to'g}: 

v/qj = R!{axv'^^ — btv’Bi) (4.43) 
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where 


a = 


1 + 


IdQ^ 


1-1 


ar 

b = a(Tf6)dQ^ldT 


(4.44) 

(4.45) 


and Qgt is the saturation specific humidity at the liquid water temperature 7). For 
now, we regard Rf as known. Substitution of equation (4.37) into equation (4.42) 
and rearrangement yields 


ti/ei = (ydi ^ 13^' 


(4.46) 


where 






(4.47) 


In order to close equation (4 37), we use equation (45) of Yamada and Mellor (1975): 



2e 


- f dQtu -y^ . ^l-TZr\ I na> Tj- 


(4.48) 


Equations (4.33), (4.34), (4.35), (4.37), and (4.46) may be regarded 
as the closure of equations (4.26)-(4.31) in that the five second-order quantities 
required to close equations (4.26)-(4.31) are expressed by these equations. 

The master length scale £ appearing in equations (4.27)-(4.29) is of the 
same form as equation (4.15), however, the asymptotic reference value is not a 
constant but rather a ratio of first to the zeroth integral moment of the turbulent 
velocity scale. The vertical integration is terminated at the top of the boundary layer 
so that isolated turbulent layers aloft cannot lead to unrealistically large values of 
£. In these layers, a constant value of 50 m is used for the length scale. 
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Details concerning the vertically nested model can be found ebewhere 
(Burk and Thompson, 1989; Thompson and Burk, 1991). Here we will only briefly 
summarize important aspects of the scheme. The computational cycle is initiated 
by taking a ‘‘dynamics time step” (St = 120 s) in which large scale advection, 
pressure gradient, and Coriolis terms are computed, yielding dynamics tendency 
terms at each point. These dynamics tendency terms provide synoptic-scale forcing 
to the prognostic set of equations which embody the parameterized physics. On 
the first model time-step, the turbulence quantities undergo a “spin-up period” of 
several iterations during which these quantities are aUowed to change while the mean 
fields and the time integration counter remain fixed. This ensures that the initial 
turbulence fields will be consistent with the regional model fields. 

The predictive equations for the parameterized physics are solved using 
implicit, tridiagonal finite difference algorithms with upper and lower boundary 
conditions inserted into the resulting matrices. A form of Gaussian elimination is 
used to achieve solution. 

As with the K-theory parameterization, no convective adjustment scheme 
is required with the second-order parameterizations. 

The three major experiments to be performed (adiabatic and inviscid, 
K-theory, and second-order closure) along with a variety of sensitivity experiments 
as discussed in Chapter III, will address the concerns highlighted at the end of 
Chapter II. 

D. EVALUATION OF FRONTOGENETIC FORCING 

The term frontogenesis is defined by Petterssen (1956, p.200) as “... a ten¬ 
dency toward formation of a discontinuity or intensification of an existing zone of 
transition...”. The intensity of frontogenetic forcing can be quantified by deter- 
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mining the rate at which the gradient oi some scalar quantity increases. These 
considerations lead to the definition of the frontogenetic function, which can be 
written as 


d\VHS\ 

dt 


(4.49) 


where S is some scalar quantity and d/dt represents a three-dimensional total deriva¬ 
tive. Following Bluestein (1986, p.181), the right hand side of (4.49) can be expanded 
and the resultant terms separated into four groups as follows; 


ds/dx r a (ds\] ds/dy r a (ds\] 
I VmS I [dx W* j J I ^hS I [ay \ dt )\ 

ds/dx (duds\ ds/dy (dvds\ 
\VHS\[dxdx) \V„S\[dydy) 

ds/dy (duds\ ds/dx (dvds\ 

\V„S\[dxdy) 

as/ax /ati;a5\ ds/dy (dwds\ 
\v„s\[dxdz) |Vh 5| \,ay azj 


[Note that Bluestein obtains additional terms in his equation (9.11) because he 
uses the threedimensional operator, V5, rather than VuS.] The first group above 
consists of gradients of the substantial derivatives of S (for the case of 5 = potential 
temperature, the first group can be identified with gradients of diabatic terms). The 
second group comprises the stretching deformation. The third group consists of the 
shearing deformation terms. The stretching deformation term quantifies the effect of 
flow normal to the isentropes in concentrating or reducing the temperature gradient 
while the shearing deformation term quantifies the effect of a horizontal gradient in 
flow parallel to the isentropes. The fourth group is referred to as the '^tilting term”. 
This term quantifies the effect of a gradient in vertical motion tilting the isentropes 
into a vertical plane. Bluestein presents further discussion and graphical depictions 
of the physical interpretation of these terms. 
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V. RESULTS 


A. ADIABATIC AND NEARLY INVISCID AT fflGH VERTICAL 
RESOLUTION 

1. Development 

In this simulation, which serves as a control for simulations discussed 
below, the parameterizations for the diabatic and turbulent processes are inactive. 
As noted in Chapter IV above, this simulation is nearly invisdd in the sense that 
horizontal diffusion, which constitutes a pseudo-viscous effect, is always present. 

The initial conditions for this case indude a linear (in z) increase in 
the u component of the wind from zero near the surface to 40 m s~^ at 200 mb 
u decreases above 200 mb. The v component consists of a sinusoidal perturbation 
with a wavelength of 3000 km and an amplitude of 1 m s'*. The geopotential field 
is in geostrophic balance with the wind field (both ti and v) and the temperature 
field is in thermal wind balance with the u wind field. Temperature decreases with 
elevation at a rate of 8 C/km. The sea surface temperature, which is held fixed 
in these simulations, is initially set equal to the temperature at the lowest model 
level. Figure 5.1 shows the sea surface temperature field. With the exception of v 
and the geopotential, all fields are constrained to be invariant in the zonal direction 
initially. Figures 5.2a and b show meridional cross sections of the initial potential 
temperature and u component, respectively (the distribution of grid points in the 
vertical can also be seen). 

Given such a baroclinically unstable initial state, rapid (^clogenesis may 
be anticipated. Figure 5.3 shows a time series of the minimum surface pressure 
for the present case (Al), the K-theory boundary layer (KTH) discussed in Section 
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5.2 bdow, and the aecond-order doaure simulation (20C) described in Section 5.3 
below. Note the rapid deepening beginning near 72 hours. By 120 h, the surface 
pressure has fallen to 965 mb. 

Figure 5.4 shows the surface pressure and near surface potential tem¬ 
perature at 72 h. Although the surface pressure has fallen only to 996 mb at this 
time, the temperature wave is well devdoped with a strong warm front and a weak, 
diffuse cold front. The thermal situation depicted in Fig. 5.4 is consistent with the 
traditional conceptual model of midlatitude cyclogenesis. Indeed, Fig. 5.4 bears a 
striking resemblance to a schematic presented by Holton (1979; Fig. 9.9, pg. 238). 
Figures 5.5a and b show the relative vorticity and divergence at approximately 850 
mb, respectively, while Fig. 5.5c shows the vertical motion field at the same level. 

Adiabatic and inviscid simulations were performed by Polovarapu and 
Peltier (1990) using a nonhydrostatic primitive equation model with periodic east- 
west boundaries and a grid structure very similar to the present one. Figure 5.5a 
resembles Fig. 13b and Fig. 5.5c closely resembles Fig. 14b of Polovarapu and 
Peltier, although Fig. 5.5 is for 72 h at 850 mb while Figs. 13b and 14b are for 
120 h at an elevation of 1 km. The relative vortidty shown in Fig. 5.5a is large 
(2.1 X lO'** s~^) within the warm frontal zone with a ‘‘tail” extending along the 
cold front. The strong warm front with large relative vorticity in the early stages 
of frontogenesis is consistent with the results of Levy (1989). The warm front lies 
within an zone of strong convergence and ascending motion (—7.8p b s~^) which 
extends towards the SW in the area ahead of the cold front. The region occupied 
by the cold air behind the front is characterized by weak divergence and descent. 

Shown in Fig. 5.6 is the near-surface potential temp^ature at hour 84. 
Substantial frontogenesis has occurred along both of the fronts and the low has 
deepened to 986 mb. Although the strong convergence is confined to the northern 
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flank of the low, the vorticity maximum extends to the southwest along the cold 
front. 

Figure 5.7 shows the near surface potential temperature for hour 96. The 
minimum pressure is 976 mb. There is a wdl defined cold front and an associated 
frontal trough. The cold front is now stronger than the warm front, although both 
fronts have experienced strong frontogenesis during the preceding 24 hours. The 
maximum in cyclonic relative vorticity has more than doubled during this period 
(not shown). Although the maximum is confined to the north, there is also strong 
relative vorticity within the cold frontal zone. An occlusion is forming near the center 
of the low and the warm air has “wrapped around” the low creating a distinctive 
pattern in the isentropes within the central portion of the low. This pattern bears a 
resemblance to the “bent back” warm front discussed by Shapiro and Keyser (1990). 
A discussion on this topic appears in the Appendix. 

At hour 108, the warm front attains its maximum strength and the 
occlusion occupies a smaller portion of the domain than it does at hour 120. Figure 
5.8 shows the near surface potential temperature. The minimum pressure at this 
time is 968 mb and the cold frontal trough is well-defined. Both of the fronts are 
quite strong (9 C (100 km)”* for the cold front and 5 C (100 km)”* for the warm 
front). 

Figure 5.9 shows the 850 mb relative vorticity field. Large values (3.6 x 
10”^ s”*) of cyclonic vorticity are confined to the inunediate vicinity of the low. The 
850 mb divergence field shows strong convergence (—0.45A^10”^ s”*) along the cold 
front. The warm front lies in a region of weak convergence. The vertical motion 
field shows similarity to the divergence field with ascending motion all along the 
cold front and with the maximum values ('v 1-2 cm s”*) to the north. There is 
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weftk descending motion in the odd air bdiind the front l e ad in g to thermally direct 
circiilation about the cold front. 


Shown in Fig. 5.10 are the surface pressure and near surface potential 
temperature at 120 h. By this time, an intense cyclone has devdoped with maximum 
near surface wind speeds of nearly 60 m s~^. The distinctive pattern in the near 
surface isentropes is more pronounced as wdl. Strong warm and cold fronts and a 
well-developed occlusion are also present. Figure 5.10b is similar to Fig. 5.10c of 
Polovarapu and Peltier (1990) for hour 144 of their /-plane/westerly jet experiment. 

2. IVontal Structure and Circulation 

A cross section of potential temperature at hour 72 is shown in Fig. 5.11. 
The cross section extends from the surface to 505 mb in the plane shown in Fig. 
5.4b. Note the relatively strong baroclinic zone at low levels in the center of the 
plane and the strongly stable surface layer to the north. The temperature of the 
surface in this region is on the order of 5 degrees wanner than the air (see Fig. 5.1), 
thus, the surface layer is strongly unstable. The surface layer becomes progressively 
less stable to the south, becoming neutral over the lowest 30 mb near the southern 
end of the cross section. The front-normal wind component (not shown) exhibits 
strong convergence near the front; it does not, however, change sign in the frontal 
zone. There is weak vertical motion near the front with ascent in the warm sector 
and descent in the cold air behind the front. 

Shown in Fig. 5.12 is the potential temperature cross section for hour 
84. The baroclinic zone is now much stronger and extends through a greater depth 
(note the region from 950 mb to 800 mb in Figs. 5.11 and 5.12). There is strong 
convergence in the front-normal component (not shown). Note that convergence 
not only strengthens the frontal temperature gradient directly but also increases 
cyclonic relative vorticity due to vortex stretching. The cyclonic circiilation im- 
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plied by increasing cyclonic vortidty acting on the along-front temperature gradient 
constitutes a frontogenetic influence for cold fronts (Gidd, 1978). Thus, increased 
convergence has both a direct and an indirect frontogenetic effect. The vertical 
circulation is slightly stronger as compared to that at hour 72. 

The potential temperature field at hour 96 is shown in Fig. 5.13 (the 
plane of the cross section is indicated in Fig. 5.7). The front, which is on the western 
boundary of the domain, has strengthened considerably. One interesting feature of 
the field is the relatively strong low levd inversion in the warm sector. Given the 
warm advection which has occurred in this re^on, high stability is not surprising. 
Nor is it surprising that the isentropes are flat as the potential temperature is 
horizontally uniform. The surprising aspect is the strength of the inversion at low 
levels. This thermal structure is the result of the combined effects of horizontal and 
vertical advection (this area is dominated by weak ascent after hour 72). Warm 
air advection results in a stable stratification. Away from the surface, vertical 
advection in the warm air results in a neutral lapse of potential temperature. As 
vertical motion decreases toward the surface, however, the positive vertical gradient 
in potential temperature remains unaltered. 

There is substantial convergence near the leading edge of the front at 
this time. A large maximum in vortidty (2/) coincides with the leading edge of the 
front but is confined to a shallow layer. 

A potential temperature cross section across the cold front at 108 hours 
is shown in Fig. 5.14 (the plane of the cross section is indicated in Fig. 5.8). Note 
the intense cross-front gradient both at the surface and above. Figures 5.15a and b 
show the components of the wind along the front and normal to it, respectively. The 
front normal wind changes sign at the front, indicating that air parcels approadi 
the front from both the cold air and the warm sector. A cross section of divergence 
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is shown in Fig. 5.16. Note that the region oi convergence extends upwards along 
the front. Given this divergence field, weak vertical would be expected (Fig. 

5.17). Note the ascent 2 cm s~^) confined to a narrow updraft immediately ahead 
of the front. The descending branch of the vertical circulation is somewhat more 
diffuse with large values some distance behind the front. 

The along-front component is strongest at the surface behind the front 
and very weak above 800 mb. This component exhibits strong cyclonic shear in the 
frontal zone. The relative vorticity is quite strong (nearly 3/ s'**) at the surface 
near the leading edge of the front (Fig. 5.18). 

3. FVontal Development 

In order to clarify the discussion of frontogeneticaJ forcing, we will first 
examine the situation at 72 h when the interpretation of the pattern is straight for¬ 
ward (refer to Section 4.4 for a discussion of the frontogenetic forcing). Figure 5.19a 
shows the total adiabatic frontogenetical forcing at hour 72. Near the warm front, 
there are two regions of strong frontogenesis separated by strong frontolysis. There 
is also frontolysis south of the warm front. The dipole in forcing with frontolysis to 
the north and frontogenesis to the south is primarily to due the titling term (Fig. 
5.19b). Naturally, the tilting term is significant only where there is ascending mo¬ 
tion. Strong frontogenesis to the east is due to the combined action of the stretching 
and shearing deformation (Figs. 5.19c and d). The region of strong frontolysis in the 
total is snutUer than the corresponding area in the tilting term due to frontogenesis 
in the stretching and shearing deformation. Weak frontogenesis along the cold front 
is due to stretching deformation in the southern portion and shearing deformation 
in the northern portion. In that portion of the front where stretching deformation is 
frontogenetic, the wind, which is not far from normal to the isentropes, is decreasing 
rapidly downstream. In that portion of the front where shearing deformation is fron- 
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togenetic, the wind turns c^clonically in the downstream direction (NW beccunes 
W becomes S). 


The total frontogenetic forcing for hour 84 (not shown) is strongly fron- 
togenetic along the warm front with very small positive values along the cold front. 
At hour 96, the distribution of the fordng changes dramatically and is strongly fron* 
togenetic along the cold front and frontolytic at the warm front. Here again, the 
frontogenesis is due to stretching along the southern portion of the front and shear¬ 
ing along the northern portion. Shearing deformation is responsible for the bulk of 
the frontolysis along the warm front. At hour 108, there is weak frontogenesis along 
the cold front with weak frontolysis along the warm front (Fig. 5.20). The tilting 
term is now predominantly frontolytic. While the stretching deformation is strongly 
frontogenetic along the southern portion of the cold front strong frontogenesis in 
the shearing deformation is confined to the northern portion of the front where the 
cyclonic vorticity is largest. There is significant cancellation between the stretching 
and shearing deformation along both the warm and cold fronts. 

A cross section of the total adiabatic frontogenetic forcing for hour 108 
(Fig. 5.21a) shows that there is strong frontogenesis both within and ahead of the 
frontal zone due primarily to stretching deformation (Fig. 5.21b). The frontogenesis 
is, in fact, exactly twice as large as at hour 96. The magnitude of the stretching 
deformation is actually larger than the total due to frontolysis in the tilting term 
superimposed on the frontal zone. The distribution of the shearing deformation is 
quite similar to the stretching deformation but it is substantially smaller in this 
plane. 

4. Summary 

The results presented in this section demonstrate that, in the absence 
of boundary layer physical processes, a growing baroclinic wave will generate very 
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strong warm and cold fronts near the surface. The numerical results are summarised 
in Table 5.1. Note that values calculated for a strong add front over the NE Pacific 
Ocean observed on 16 Novend>er 1980 as doctimented by Fleagle et (1988) are 
included in the last row for comparison. The rdative vortidty and convergence 
within the frontal zone are large, as is the horizontal potential temp^ature gradioit 
across the front. The results are in qualitative agreement with other published 
studies of adiabatic and invisdd frontogenesis. 

These simulations, however, are unrealistic in several ways. For example, 
the wind speeds near the center of the low at the surface are too large (~ 60 m s~^). 
The vertical motion and the circulation about the front are extremely weak; the 
ascent ahead of the front is equivalent to 1~2 cm s~*, nearly an order of magnitude 
below observed vertical motion in the vicinity of cold fronts. In contrast to previous 
kinematic studies of frontogenesis, tilting frontogenesis is of little importance in the 
total frontogenetic forcing. 

In the next section, boundary layer processes are included using the K- 
theory parameterization discussed in Chapter IV above. Addition of surface drag 
should both reduce the low-level wind speeds and strengthen frictional convergence 
which will increase upward motion. Turbulent mixing of heat and momentum will 
also be active. Thus, the results should be more representative of observed fronts. 

B. K-THEORY BOUNDARY LAYER 
1. Development 

In this experiment, the initial conditions are identical to those in the AI 
case (Section 5.1 above). The K-theory parameterization is activated at hour 72. 
Thus, the evolution of the baroclinic wave for the first three days is identical to that 
discussed in the AI case. A time series of surface pressure is shown in Fig. 5.22 


50 








for the present case (KTH), the simulation with surface sensible heat flux excluded 
(NO FLUX) described in Section 5.2.4 below, and for the simulation with neutral 
eddy coefficients (NO STABILITY) described in Section 5.2.5 below. Note that less 
deepening occurs after hour 72; the mininmm surface pressure decreases 10 mb in 
the following 48 hours (in the AI case, the pressure decreased by 30 mb over this 
period). Time aeries of minimum surface pressure for the present case and the AI 
case are shown together in Fig. 5.3. 

The near-surface potential temperature 12 hours after the K-theory pa- 
rameteri 2 sation is activated (84 hours) is shown in Fig. 5.23. The minimum pressure 
at this time is 990 mb and the cold and warm fronts are relatively weak. 

Figure 5.24 shows the near surface potential temperature at hour 96. At 
this time, there is a well defined frontal trough near the cold front. Both fronts 
are relatively weak, although the cold front is slightly stronger than the warm front 
(2 K (100 km)”^ for the cold front; 1.7 K (100 km)”* for the warm front). The 
divergence field at 850 mb (not shown) indicates weak convergence along the cold 
front. Weak divergence dominates the high pressure region. The vertical motion 
field exhibits a narrow region of weak ascent extending along both fronts with weak 
descent bdiind the cold front. 

Figure 5.25 shows the near surface potential temperature for hour 108. 
Under the influence of continuing frontogenesis, the warm and cold fronts have 
intensified. The frontal trough is also well developed. The 850 mb relative vorticity 
is cyclonic along both fronts, although it is smaller along the warm front than it is 
along the cold front (Fig. 5.26). 

At hour 120, the minimum central pressure has fallen to 984 mb and 
the cold front has strengthened considerably (Fig. 5.27). The warm front has 
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experienced frontolyais. The maxiinuin near surface wind speed is 23.6 m s~* (as 
compared to nearly 60 m s~^ in the AI case at 120 h). 

The baroclinic wave amplitude is smaller in this case and the potential 
temperature gradients associated with the warm and cold fronts are dramatically 
different from the AI case. There are several reasons for this. The wind speed is 
lower due to surface drag. Temperature gradients and, therefore, available potential 
energy, are reduced due to upward surface heat flux in the cold air. Frontogenesis 
is reduced since the geostrophic horizontal deformation is weaker. Branscome et al. 
(1989) investigated the effect of surface heat and momentum fluxes on the nonlinear 
development of baroclinic waves. They found that nonlinear surface drag was most 
important in reducing wave amplitude relative to the inviscid case. This resulted 
from enhzinced dissipation. Reduced wind speeds were also responsible for most 
of the reduction in horizontal heat transport at low levels. An interesting result 
is that the inclusion of surface heat flux causes the maximum eddy temperature 
variances to occur at 800 mb rather than at the surface as in the inviscid case. 
The results of Branscome et al. are consistent with the present study. The wave 
amplitude is substantially reduced relative to the AI case and, as will be shown 
below, the maximum frontal gradients are displaced vertically and occur at 700-800 
mb. Simulations discussed below in which surface heat and momentum fluxes are 
removed are also consistent with the results of Branscome et al.. 

2. FVontal Structure and Circulation 

A cross section of potential temperature at hour 84 is shown in Fig. 5.28. 
While there is some evidence of frontal structure in the temperature field, the front 
at this time is quite weak. Comparison of this figure to the corresponding one for 
the AI simulation (Fig. 5.12) clearly shows the influence of surface heat flux; the 
strongly unstable stratification in the cold air in the AI simulation is replaced here 
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by a 150 mb deep mixed layer. Strong surface fluxes in the unstable air behind 
the front promote entrainment and deepening of the boundary layer in the present 
case. 

A cross section of potential temperature for hour 96 is shown in Fig. 

* 5.29. By this time, the cross-front temperature gradient is relatively strong and the 

ageostrophic circulation about the front has become fairly vigorous. The region of 
positive relative vortidty corresponds closely to the location of the largest horizontal 
temperature gradient. The relatively low stability in the warm sector is the result 
of weak ascent in this region beginning between hours 84 and 96. The frontal zone 
in this case is dramatically different from the AI simulation at this time (Fig. 5.13), 
although some similarities are apparent in the warm sector. In this case, the front 
is weak in comparison to the AI case and there is a deep mixed layer in the cold 
air. Reduction in the potential temperature gradient across the front in the present 
case is the result of strong upward heat flux in the cold air. 

A cross section of potential temperature at hour 108 is shown in Fig. 
5.30. Comparison to F'g. 5.29 shows that the front has moved approximately 500 
km over the preceding 12 hours. Note the well-mixed layers both ahead of and 
behind the cold front extending to 800 mb. In the frontal zone, there is a weak 
inversion extending from 800 mb to 700 mb with a much stronger inversion above. 
In the cold air behind the front (north of the 295 K isentrope), there is a much more 
well-defined inversion at 700 mb. In the strong inversion above the frontal zone, 
the horizontal potential temperature gradient is twice as large as the near-surface 
gradient. The potential temperature distribution in the AI case shown in Fig. 5.14 is 
radically different from the present case. The frontal potential temperature gradient 
is less than half as large as in the AI simulation and the mixed layer is ~ 200 mb 
deeper. Clearly, the surface fluxes of heat and momentum and turbulent diffusion 
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have had a large impact on frontal structure and circulation. Figure 5.31 shows 
cross sections of the along-front and front-normal wind components. Note the low 
ievel-jet in the cold air and strong horizontal shear in the along front component 
and strong convergence in the front-normal c.>mponent near the front. The region in 
which convergence is large extends over a layer approximately 150 mb deep. There 
is very little vertical shear over the lowest 150 mb, indicating that momentum is 
well mixed in this layer. This is in contrast to the AI case in which the maximum 
along front wind is at the surface and the front-normal wind exhibits substantial 
vertical shear at low levels (Figs. 5.15a and b). The convergence (Fig. 5.32) near the 
front extends through this layer and back along the frontal zone. There is relatively 
strong divergence above the surface convergence. The strong convergence results 
in relatively strong upward motion (Fig. 5.33). There is an ascending jet at the 
leading edge of the front with strong descent at upper levels in the cold air. The 
strong cyclonic shear in the along-front component associated with the low-level jet 
in the cold air leads to a maximum in cyclonic relative vorticity in the frontal zone 
(Fig. 5.34). 

An important result is that the convergence is only about half the value 
in the AI case (0.4 x 10“^ s~^ vs. 0.7 x 10“^ s~* for the AI case) but that, near 
the leading edge, the convergence extends through a deep layer. This is evident in 
comparing Figs. 5.16 and 5.32. This convergence pattern gives rise to an updraft 
more than twice as strong as in the AI case (4.4 cm sec”^ vs. 2.0 cm sec“^ for 
the AI case). The ascending jet is also smaller in both its horizontal and vertical 
extent while the descending branch of the ageostrophic circulation is also more 
compact. Thus, the turbulent mixing in concert with the surface fluxes organizes 
the convergence in such a way as to produce a much more compact and vigorous 
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transverse circulation. This is primarily the result of mixing of momentum. This 
will be discussed at length in the following chapter. 

It is of interest that the frontal zones depicted in the potential temper¬ 
ature fields displayed above in Figs. 5.29 and 5.30 exhibit no tilt in the vertical 
over the lowest 100-200 mb. This is in contrast to the adiabatic and nearly invisdd 
simulations and many observational studies (e.g., Sanders, 1955; Ogura and Portis, 
1982) in which the frontal zone has a pronounced vertical tilt. The lack of tilt results 
from the deep well-mixed boundary layers on both sides of the front. This is to be 
expected in the cold air, however, one might anticipate higher stability in the warm 
sector. Nevertheless, the simvdated frontal structure is consistent with the structure 
of observed marine cold fronts (Fleagle et al., 1988). This will be discussed in detail 
below in the following chapter. 

3. Frontal Development 

The total adiabatic frontogenetic forcing at hour 84 at 850 mb shows 
strong frontogenesis and frontolysis along the northern portion of the warm front. 
This is primarily due to tilting with large frontogenetic contributions from shearing 
and stretching deformation. At hour 96, there is strong frontolysis along the cold 
front and behind the northern portion of the warm front. This is almost entirely 
due to tilting with stretching and shearing deformation contributing only along the 
northern portion of the warm front. The total adiabatic frontogenetic forcing at 
850 mb at hour 108, shown in Fig. 5.35, is predominantly frontolytic along both 
the warm and cold fronts. There is weak frontogenesis ahead of both fronts. This 
is largely due to tilting, although the frontolysis along the cold front is partially 
canceled by stretching deformation to the south and shearing deformation to the 
north. Neither of the horizontal deformation terms contributes significantly along 
the warm front. 


55 




Figure 5.36 shows a cross section of the total adiabatic frontogenetic 
forcing at 84 hours. Note the dipole structure of the forcing in the layer centered at 
approximately 800 mb. This pattern is due almost entirely to the tilting term (not 
shown); the shearing and stretching deformation are both an order of magnitude 
smaller than the tilting term. The distribution of the tilting term is a consequence 
of the vertical motion field (not shown) whidx consists of a weak downdraft in the 
cold air behind the front. The horizontal gradients in vertical motion associated with 
the downdraft result in frontogenesis well behind the front with weak frontolysis in 
the frontal zone extending into the warm sector. 

The frontogenetic forcing at 96 hours is shown in Fig. 5.37. The fron¬ 
tolysis at the leading edge of the front is due to tilting. At this time, the vertical 
motion field consists of a relatively strong updraft with a large graidient in vertical 
motion on the cold side and a broad, diffuse region of upward motion on the warm 
side. This results in strong frontolysis within the frontal zone and weak frontogen¬ 
esis in the warm sector. The total frontogenetic forcing shows a dipole with strong 
frontogenesis on the cold side of the tilting frontolysis. This is due to the shearing 
deformation, which is strongly frontogenetic both at the surface at the leading edge 
and above and behind the leading edge. 

Figure 5.38a shows a cross section of the total adiabatic frontogenetic 
forcing at hour 108. Note the dipole structure with strong frontolysis within the 
frontal zone and strong frontogenesis immediately ahead of the front centered at 
850 mb. This is almost entirely due to the tilting term (Fig. 5.38b). It is an 
order of magnitude larger than the stretching and shearing deformation and an 
order of magnitude larger than the tilting term in the AI case. Near the surface, 
where the tilting term is small (it drops to 50% of its maximum value at 950 mb and 
decreases rapidly toward the surface), the stretching and shearing deformation play a 
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more important role. Both terms have a maxima near the surface and the shearing 
deformation extends upward along the frontal surface. The maximum values are 
near the surface and above the inversion where the horizontal temperature gradient 
is largest. 

At this time, the updraft ahead of the front is much more compact and 
symmetric than at hour 96 (Fig. 5.33). The dipole structure apparent in both 
the tilting frontogenesis and the total fronto^nesis results from this distribution of 
vertical motion. The large gradients in vertical motion result in large magnitude 
in the tilting term. The largest gradient in vertical motion is within the frontal 
zone where the vertical motion changes sign. The tilting term is strongly frontolytic 
within the frontal zone. Ahead of the front, on the warm side of the vertical jet, 
the reversal in direction of the gradient in vertical motion results in a frontogenetic 
contribution from the tilting term. 

At hour 108, the frontal propagation speed is larger than the front- 
normal wind speed, which decreases to zero in the frontal zone and then changes 
sign in the warm sector. This indicates that, relative to the front, air parceb in the 
cold air are moving away from the front while air parceb in the warm sector are 
moving towards it. Thus, air parceb in the cold ak are moving away from the region 
in which the temperature gradient is large while air parceb in the warm sector move 
towards it. This is consistent with a dipole in tilting frontogenesis from a Lagrangian 
point of view. 

Figure 5.38c shows the diabatic frontogenetic forcing. It is strongly fron¬ 
tolytic at the leading edge of the frontal zone and frontogenetic at the trailing edge. 
The magnitude of the frontolysb exceeds the frontogenesis and the region occupied 
by frontolysis is larger. If the diabatic term were the only term operating, the front 
would be destroyed. 
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The diab&tic tenn involves at least two ui 4 >ortant sources of frontolysis 
in this case. One is the frontolytic effect of vertical mixing. In the cold air behind 
the front, there is a deep mixed layer capped by a strong inversion. Entrainment 
of warm air into the cold mixed layer has the effect of reducing the temperature 
gradient across the front. Similarly, heat flux divergence associated with upward 
surface heat flux in this region will reduce the cross-frront temperature gradient. 
The distribution of surface sensible heat flux (W m~^ ) for hour 108 is shown in 
Fig. 5.39. Note that there is upward surface heat flux in the cold air behind the 
cold front. The effects of surface heat flux and entrainment are not independent in 
a well mixed BL. If the effect of the surface sensible heat flux were removed from 
the simulation, the frontolytic influences of both surface heat flux and entrainment 
should be reduced. Due to the highly nonlinear interaction between the surface 
heat flux and the parameterized turbulent mixing, however, differences between 
simulations with and without surface heat flux will involve substantial modifications 
in the turbulence as weU. 

Turbulent mixing is dependent upon wind shear and static stability as 
well as surface heat flux. In the next subsection, results from a simulation in which 
surface heat flux is eliminated are discussed. In Section 5.2.5, the effect of static 
stability is removed from the eddy coeffidoits. 

4. Insulated Surface Simulation 

The near surface potential temperature field at hour 84 of the simulation 
with the surface heat flux removed is shown in Fig. 5.40. The warm and cold fronts 
are much stronger in this case. Over the next 12 hours, the minimiim pressure falls 
4 mb and both fronts experience significant frontogenesis (Fig. 5.41). Frontogenesis 
continues over the following 12 hours. Figure 5.42 shows the near surface potential 
temperature for hour 108. While the pressure distribution is nearly the same as in 


58 











the full physics simulation, the potential temperature fidd is quite dififerent. The 
warm and cold fronts are both much stronger. The warm ‘‘tongue" ahead of the 
cold front is much more pronounced. 

The fact that the surface pressure distribution is similar to that for the 
simulation with heat flux included indicates that surface heat flux has little impact 
on baroclinic wave growth. This result is consistent with the study of Branscome et 
al. (1989) discussed above. A time series of nunimum surface pressure is shown in 
Fig. 5.22 with the time series for the full pl^cs simulation for comparison. 

Shown in Fig. 5.43 are cross sections of the potential temperature for 
hours 84,96, and 108. Note that the front in this case has a very small vertical scale 
as compared to the full physics simulation. The front perturbs only the lowest 30 mb 
of the atmosphere. Since the influence of the front extends over such a limited depth, 
there is very little response to it in terms of vertical motion. The upward motion at 
hour 108 amounts to only about 1 cm s*’^ There is frontolysis in the tilting term 
but it is a result of large horizontal gradients in ut which are in turn the result of 
the rather disorganized structure of the ut field; this frontolysis plays no significant 
role in the frontal kinematics. The horizontal deformation is frontogenetic, but it is 
confined to a shallow depth. Thus there is no dipole structure in the frontogenetic 
forcing; this case is, in this respect, similar to the AI case. 

The present case is similar to the AI case in several respects, although 
there is a 30 mb deep mixed layer in the present case which is absent in the AI 
case. The total change in potential temperature across the front, for example, is 
nearly the same as in the AI case. The width of the frontal zone, however, is twice 
as large in the present case. In both cases, the vertical motion is very weak. 

The temperature gradient across the front is much larger in this case 
than in the full physics simulation. The enhanced cross front temperature gradient 
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reflects the frontolytic effect of suiihoe fluxes in this case, while the shallow mixed 
layer implies that the entraininmit has been reduced by excluding the heat flux. The 
mixed layer is maintained by vertical wind shear in the front-normal component. 

5. Neutral Eddy CoeflScient Simulation 

The form of the eddy coefficients in the model is given by 




dV 


dz 


F{IU) 


(5.1) 


where L = kz/{l -f- kz/X) and where A s 150 for and 350 for Km (for definitions 
and a more complete description, see Section 4.3.3). Dependence on stability can 
be eliminated by removing F(Ri); the neutral value of the eddy coefficient is thus 





av 

dz 


(5.2) 


Figures 5.44a and b show the normalized eddy momentum and heat coefficients, 
respectively, as a function of Richardson Number for a constant wind shear of 10 m 
s~^ km**. The line having zero slope is the neutral value. Under stable conditions 
{Ri > 0), the neutral value is, not surprisingly, a considerable overestimate of the 
eddy coefficient, while, for unstable conditions {Ri < 0), the neutral value is a 
somewhat better estimate of the stability dependent value, particularly for the eddy 
heat coefficient. 

Owing to the stabilit b \ acteristic of the free atmosphere away from 
the surface, the formulation (5.1) gives small values for the eddy coefficient above 
the BL. In a simulation using (5.2), large eddy coefficients at high elevations lead to 
unphysical results. Therefore, (5.2) was modified to include a parameter dependent 
on <r: 


Kk,m = 


dV 


dz 


(5.3) 
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where 


0 = 1 1 < o < 0.85 

o s 0.1 0 < O' < 0.85 

This formulation, though somewhat arbitrary, gives satisfactory results. 

The near surface potential temperature distribution for hour 108 is shown 
in Fig. 5.45. While the surface pressure distribution is similar to the full physics 
case (see Fig. 5.22), the potential temperature distribution is quite different (see 
Fig. 5.25). Ahead of the cold front, the “tongue" of warm air extends much further 
to the NW. This takes place in a repon whidi is stable (see the surface heat flux 
distribution in Fig. 5.39 for the full physics simulation) and is the result of enhanced 
turbulent mixing (recall that the neutral eddy coeffident overestimates the stability 
dependent eddy coefficient in stable layers). Surface fluxes are stiU active in this 
case, cooling the near surface layer. However, increased mixing of warm air toward 
the surface reduces heat flux divergence. A similar feature appears in Fig. 5.42 for 
the no heat flux case. 

Shown in Fig. 5.46 is the potential temperature cross sections in the 
vidnity of the cold front at hoiu: 108. The cross sections at hours 84 and 96 (not 
shown) are similar to those in the full physics simulation at corresponding times. 
Although the strength of the front at hour 108 is unchanged, the well defined in¬ 
version and mixed layer evident in Fig. 5.30 are absent in this simulation. This 
is to be expected given the reduced eddy coefficients in the unstable region in the 
vidnity of the front. The extremely unstable lapse near the surface is also a result 
of reduced mixing. The along-front and front-normal wind components (not shown) 
exhibit larger vertical shear above the BL and more divergence in the cold air with 
a slight increase in convergence near the front. Momentum is well mixed only to 850 
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mb, which is not surprising given the specification of o in equation (5.3). The alight 
increase in convergence over a slightly reduced depth results in ascmding motion 
of the same magnitude as in the full physics simulation. The region of descent is 
somewhat larger due to enhanced divergmce in the cold air. 

6. Summary 

Use of the K-theory parametoization resulted in a more realistic firontal 
structure than in the AI case. The numerical results are summarized in Table 5.1. 
Note that the values calculated for a strong cold front over the NE Pacific Ocean 
observed on 16 November 1980 as documented by Fleagle et al. (1988) are included 
in the last row for comparison. The vertical circulation was much stronger, in spite 
of the fact that the convergence was not as strong. Thb was due to the action of 
parameterized BL processes. 

Simulations in which various physical processes were excluded indicated 
that surface heat flux was most important in developing frontal structure. Although 
the direct eff^t of surface heat flux was frontolytic, the indirect effect through 
modulation of BL structure was critical in producing realistic frontal features. A 
simulation was performed in which surface drag was removed. Removal of surface 
drag resulted in stronger cyclonic circulation and higher frontal propagation speed 
which enhanced convergence and vertical motion ahead of the front. When the effect 
of static stability was removed from the eddy coefficients, details of BL and frontal 
structure were lost but the gross features were unchanged. 

Both the convergence and vorticity in the frontal zone were smaller than 
in the AI case. In a simulation in which both surface heat and momentum fluxes were 
removed, convergence and vorticity were nearly doubled and were then comparable 
in magnitude to the AI case. 
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In the next section, boundary layer physical processes are included using 
the second-order closure turbulence parameterization discussed in Chapter IV above. 
While the results described in the previous section demonstrated the impact of a 
variety of boundary layer processes on frontogenesis through comparison to the 
adiabatic and invisdd case, the more realistic treatment of turbulence afforded by 
the second-order closure scheme will serve to highlight the specific processes most 
important in maritime atmospheric frontogenesis. 

C. SECOND-ORDER CLOSURE BOUNDARY LAYER 
1. Development 

In this simulation, the second order closure physics is activated at hour 
72. Very little deepening occurs after this time. The minimum pressure at hour 
120 is identical to that in the K-theory simulation. Figure 5.47 shows time series of 
minimum surface pressure for this case (20C), for a simulation in which the surface 
heat flux is removed (20C NO FLUX) described below in Section 5.3.4, and for a 
simulation in which the surface drag is removed (20C NO DRAG) described below 
in Section 5.3.5. Figure 5.3 shows the time series of minimum surface pressure for the 
present case, the K-theory simulation (KTH), and the AI simulation (AI) together 
for comparison. Figure 5.48 shows the near surface potential temperature 12 hours 
after the physical parameterizations are activated. The low has deepened slightly at 
this time and the center and frontal trough are better defined than in the K-theory 
simulation. Temperature gradients at both the warm and cold fronts are also larger. 
By hour 96, the cold frontal trough has become still more pronounced and the 
minimum pressure is lower (992 mb for the K-theory vs. 988 mb for the second- 
order closure). The warm and cold front are also stronger (Fig. 5.49). Figure 5.50 
shows the near-surface potential temperature at hour 108. Note the strength of the 
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fronts. The surface pressure frontal trough (not shown) has sharpened considerably 
during the preceding 12 hours as compared to the K*theory simulation. The 850 
mb relative vortidty at hour 108 is shown in Figure 5.51. The strongest cyclonic 
vorticity is associated with the low. The vorticity along the <x>ld front is larger than 
in the K-theory simulation. Figure 5.52 shows the surface pressure and near surface 
potential temperature for hour 120. The low is somewhat deeper at this time. Note 
the sharp cold frontal trough and well-defined warm frontal trough. 

The strength of the warm front and sharpness of the frontal troughs are 
unique to this simulation. These features are the result of strong convergence at the 
front. 

2. FVontal Structure and Circulation 

Figure 5.53 shows a cross section of potential temperature at hour 84. 
There is a well-defined frontal zone at this time which is nearly vertical below 930 
mb. Above 900 mb, the frontal slope is not unlike the slope of typical midlatitude 
cold fronts. The strength of the potential temperature gradient across the front 
is approximately 2 C (100 km)~*, or roughly twice as strong as in the K-theory 
simulation. The boundary layer in the cold air is not as well mixed as in the K- 
theory simulation. 

During the 12 hours from 84 to 96, the front undergoes significant fron- 
togenesis (Fig. 5.54); the cross front potential temperature gradient reaches 4 C 
(100 km)~^. Frontogenesis also occurs in vorticity (Fig. 5.55). The frontal zone is 
now vertical to ^ 890 mb and has a larger slope at higher levels as well. The basic 
structural features of the front are similar to those for the K-theory simulation, al¬ 
though the front is stronger and vertical over a smaller depth. Strong convergence 
is apparent throughout most of the cross section. The most prominent feature of 
the v^tical motion field, shown in Fig. 5.56, is the ascending jet above the surface 
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front with a maximum of 13 ^ b s~^ ('^ 12.5 cm 8~^). This value is approximately 
twice as large as values reported by Fleagle et al. (1988) over the NE Pacific ocean 
and Ogura and Portis (1982) over the south-central U.S., but only half as large as 
reported by Sanders (1955) over the south-central U.S.. Most importantly, it is also 
nearly three times as large as in the K-theory simulation. 

There are, in fact, three relative maxima in ascending motion in the 
warm sector shown in Fig. 5.56. It should be noted in passing that the three 
regions of ascent ahead of the front are somewhat suggestive of frontal rainbands. 
The broad region of ascent in the warm sector at this time in the K-theory simulation 
is a poorly resolved depiction of a similar feature. In view of the relatively coarse 
horizontal resolution used in these simulations and the lack of moist physics, the 
significance of this result is unclear. 

The potential temperature cross section for hour 108 is shown in Fig. 
5.57. The frontal structure at this time is very similar to that at 96 hours, although 
the front has moved toward the south at an average speed of 12 m s~^ (note that the 
plane of the cross section is the same as at hour 96). This structure is similar to that 
observed in the Storm Transfer and Response Experiment (STREX) documented by 
Fleagle et al. (1988). In particular, the cross section through the strong maritime 
cold front of 16 November 1980 (page 55, Fig. 4 of Fleagle et al.) shows many 
of the same features seen in Fig. 5.57 (e.g., the deep mixed layer and the strong 
horizontal gradient across the front above the boundary layer). The results of the 
K-theory simulation do not compare as favorably to the STREX results in that 
the depth of the mixed layer in the second-order closure simulation is closer to the 
observations. Figure 5.58 shows the along-front and front-normal wind components 
at hour 108. Strong convergence and cyclonic vorticity are apparent in the wind 
field. Note the sign reversal in the front-normal wind within the frontal zone and 
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the strcmg cyclonic shear in the along-front wind on the cold side of the low>level 
jet in the cold air. Figure 5.59 shows a cross section of the relative vorticity and 
Figure 5.60 shows the divergence for this time. The vertical motion field is shown 
in Fig. 5.61. Note that the maximum ascent is nearly as large as at hour 96 but 
• that the region of significant upward motion is much smaller (compare the 30 mb 
hr~^ contour in figures 5.56 and 5.61). 

The near-surface winds in this simulation are much stronger than in the 
K-theory simulation, indicating that the momentum from higher levels in the bound¬ 
ary layer is being mixed to the surface, thereby strengthening the wind speed and 
thus the horizontal shear, vorticity, convergence, and vertical motion in comparison 
to that simulation. 

In summary, the comparison of the present results to those for the K- 
theory simulation has shown that, at hour 84 (Figs. 5.28 and 5.53) the mixed 
layer in the K-theory simulation is much more well defined, although the frontal 
gradient is weaker than in the second-order closure simulation. The BL in the 
later simulation is actually slightly stable at this time. In terms of gross structural 
features, however, the two fronts are similar. Comparison of Figs. 5.29 with 5.54 
reveals that, at this time (hour 96), the size of the temperature gradient across the 
front is the most significant difference between the two simulations. The convergence 
is stronger in the second-order closure simulation and the upward vertical motion 
is nearly three times as large. Cross sections of potential temperature from the two 
simulations at hour 108 are shown in Figs. 5.30 and 5.57. The fronts at this time 
are radically different. The cross-front temperature gradient is more than twice as 
large in the second-order closure simulation. The convergence and vertical motion 
are also stronger in this simulation due to more vigorous mixing of momentum. The 
cyclonic shear across the front in the along-front wind component is twice as large 
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in thb simulation, giving rise to much larger cyclonic vortidty. The post-firontal 
mixed layer is 150 mb deep in the present simulation and nearly 300 mb deep in the 
K-theory simulation. 

3. fVontsd Development 

The total adiabatic frontogenetic forcing for hour 108 is primarily £ron- 
tolytic along the warm front and there is frontogenesis along the cold front with 
strong frontolysis behind it (Fig. 5.62). The bulk of the frontogenesis is due to 
tilting and nearly all of the frontolysis results from tilting. The stretching deforma¬ 
tion is frontolytic along the warm front and frontogenetic along the cold front. The 
shearing deformation is 0 along the warm front and frontogenetic along the cold 
front. The distribution of frontogenetic forcing is not as well organized as in the 
K-theory simulation (Fig. 5.35). 

The frontogenetic forcing at hours 84 and 96 is shown in Figs. 5.63 and 
5.64. At these times, there is strong frontogenesis within and behind the frontal 
zone with weak and disorganized frontolysis ahead of the front. The frontogenetic 
forcing at hour 108 is shown in Fig. 5.65a. At this time, the vertical motion 
field shows a compact, well defined vertical circulation with a single ascending jet 
immediately ahead of the front and descent behind it. The tilting term shows a 
dipole signature (Fig. 5.65b). The dipole structure is n. t as apparent in the total 
frontogenetic forcing because both shearing and stretching deformation cancel the 
frontolytic effect of the tilting term. A second region of frontogenesis behind the 
front is due to shearing deformation. The total fordng is frontogenetic although 
the net effect of tilting is frontolytic. In the K-theory simulation at this time, 
while the magnitude of the total forcing is much smaller, the dipole signature is 
evident (Fig. 5.38) since the shearing and stretching deformation are both an order 
of magnitude less than the tilting term. The greater significance of stretching and 
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shearing deformation in the present simul^ion is a reflection of stronger convergence 
and vorticity. 

Viewed in this context, it might be concluded that the overall effect of 
the vertical circulation is to destroy the front (since the net effect of the tilting term 
is frontolytic in both the K-theory and second-order closure simulations). While it is 
true that the vertical circulation acts to restore geostrophic balance, consideration of 
the effect of the concentrated, strong updraft ahead of the front leads to a different 
conclusion. The updraft enhances frontogenesis both by the frontogenetic effect of 
tilting ahead of the front and because strong convergence at the base of the updraft 
increases the temperature gradient through stretching deformation and concentrates 
cyclonic vorticity at the front. Vorticity is also generated through convergence and 
by the tilting effect on vorticity. Since the along-front southerly flow in the warm 
sector decreases with elevation, horizontal vorticity is generated in the front-normal 
plane. With the gradient in vertical motion on the warm side of the vertical jet, 
cyclonic vorticity is tilted into the vertical. With the front-relative front-normal wind 
directed toward the front, this vorticity will be concentrated at the front. Thus, the 
strong updraft is associated with frontogenesis both in temperature (through tilting 
and stretching deformation) and in vorticity. 

Given that both the ascending jet and the horizontal gradient in vertical 
motion are much stronger in this case than in the K-theory simulation, the impacts 
of the updraft in terms of tilting frontogenesis, stretching deformation, concentra¬ 
tion of vorticity, and vorticity generation through tilting are more apparent in this 
simulation. 

The diabatic frontogenetic forcing is strongly frontolytic throughout the 
period. The maximum frontolysis is always located at the leading edge of the front 
and occurs at hour 96. Diabatic frontogenesis occurs within the frontal zone near 
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850 mb at hours 96 and 108. This location is near the top of the mixed layer whore 
entrainment of high potential temperature air leads to warming. It appears that the 
inversion is stronger behind the front than it is in the warm sector. This distribution 
of inversion strength indicates that there is greater warming due to entrainment in 
the warm sector, which constitutes a frontogenetic effect. This effect was not present 
in the K-theory simulation. 

For additional discussion on interpretation of frontogenetic forcing, see 
the frontal development section for the K-theory simulation above. 

4. Insulated Surface Simulation 

Shown in Fig. 5.66 is the surface sensible heat flux for hour 108. Compar¬ 
ison of the surface sensible heat flux distributions for the K-theory and second-order 
closure simulations provides a striking contrast. Examination of the heat flux distri¬ 
bution for the K-theory simulation (Fig. 5.39) shows a very weak horizontal gradient 
in the vicinity of the front. The heat flux reaches 100 W m“* approximately 400 
km behind the front in the cold air. In the second order closure simulation, there is 
an extremely strong gradient across the cold front with a maximum of 150 W m~* 
adjacent to the front on the cold side. The difference across the front is in excess of 
120 W m~^. There is a secondary maximum of 180 W m“* on the southwest flank of 
the low where the air-sea temperature difference is largest in both simulations. The 
gradient across the warm front is negligible in the K-theory case and 80 W m~^ for 
the present simulation. Given such a large disparity between the two simulations 
in surface heat flux, it is of interest to compare K-theory and second-order closure 
simulations with the heat flux removed. 

The near surface potential temperature for hour 108 of the insulated 
surface simulation is shown in Fig. 5.69. Figures 5.67 and 5.68 show the potential 
temperature fields for hours 84 and 96, respectively. While the surface pressure 
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distribution is similar to that for the full physics case (see Fig. 5.47), the distribution 
of the near surface temperature is quite different. The temperature is lower, the 
warm air has “wrapped around” the low as in the AI case, and the cold front is 
stronger than in the full physics case. The warm front, however, is not as strong. 

' Shown in Fig. 5.70a is a cross section of potential temperature at hour 

84 in the vicinity of the cold front. The front at this time is weak and the structure 
is very s imil ar to that for the full physics simulation (Fig. 5.53). The wind and 
vertical motion distributions are also similar, although the maximum upward motion 
is displaced upward by 70 mb. At hour 96, the warm sector is much cooler and 
the mixed layer behind the front is not as deep (Fig. 5.70b). The wind and vertical 
motion fields are similar to those in the full physics case. 

Figure 5.70c shows a cross section of potential temperature at hour 108 
of the simulation with surface heat flux excluded. Note that the depth of the mixed 
layer is much less in this case and that the cold air behind the front is much colder, 
although the strength of the front is nearly the same as in the previous simulation. 
There is much more vertical wind shear in the front-normal component above ~ 950 
mb in this case. The maximum convergence is displaced vertically to near 850 mb 
and is to the rear of the surface front. The maximum upward motion, which is cen¬ 
tered near 750 mb, is due to this elevated maximum and is thus completely unrelated 
to the surface convergence. This indicates that frontogenetic forcing associated with 
tilting will be insignificant at low levels, as in both the K-theory simulation with 
heat flux removed and the AI simulation. 

The intriguing aspect of this case is that, although it is slightly stronger, 
the front has not been substantiaUy strengthened by the removal of the surface heat 
flux. Both the cold air behind the frx>nt and the warm sector are 2-3 C colder than 
in the full physics simulation. Thus, the temperature difference between the cold 
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air and the warm sector is nearly tL.c same. This is also the case for the K-theory 
simulation with heat flux removed. In that simulation, however, the horizontal 
extent of the frontal zone decreases substantially, giving rise to a much stronger front 
(in fact, the strength of the front in the two simulations with heat flux removed is the 
* same). Apparently, removal of surface heat flux allows the cold air to flow unmodified 
to the rear of the frontal zone, just as in the second-order closure simulation. The 
primary difference between the two simulations is in the depth of the mixed layer. 
In the second order closure simulation, the mixed layer is approximately half as deep 
as in the full physics simulation ('>' 80 mb vs. 150 mb). In the K-theory simulation, 
the mixed layer depth is an order of magnitude less deep (~ 30 mb vs. 300 mb for 
the full physics case). 

In view of the much stronger heat flux in the second-order closure sim¬ 
ulation, its removal might be expected to have a larger effect than a removal in the 
K-theory simulation. Recalling, however, the similarities between the no heat flux 
K-theory simulation and the AI simulation leads to the observation that removing 
heat flux from the K-theory simulation is tantamount to removing all boundary layer 
processes from the K-theory parameterization. This is because surface heat flux is 
the most important process in deepening of the boundary layer. In the second-order 
closure simulation, other physical processes play a significant role as well. 

In the next subsection, the effect of surface drag is examined. 

5. Free Slip Surface Simulation 

Figure 5.73 shows the near surface potential temperature for hour 108 of 
the free-slip surface simulation. Figures 5.71 and 5.72 show the potential tempera¬ 
ture fields at hours 84 and 96, respectively. The near surface potential temperature 
is quite similar to the full physics simulation (Fig. 5.50). The low is significantly 
stronger in this case, as it was in the K-theory simulation with no surface drag. A 
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time series of minimum surface pressure is shown in Fig. 5.47 (the time series for 
the insulated surface and full physics simulations are also shown for comparison). 
The frontal trough is not as well-defined in this case, however. 

Examination of cross sections of potential temperature for hours 84, 96, 
and 108 shows that only at hour 108 are there significant differences between these 
cross section and those for the full physics case. The potential temperature cross 
sections for the present case and the full physics case at hour 108 are shown in 
Figs. 5.74 and 5.57). In this case, the mixed layer depths are slightly bwer on 
both sides of the front. The wind components (not shown) are substantially higher, 
as might be ocpected, but the convergence is only slightly stronger (Fig. 5.75). 
The vertical motion is somewhat smaller at hour 84 but stronger at hour 96 and 
substantially stronger at hour 108 (Fig. 5.76). The maximum upward motion is 
at roughly the same level as in the full physics simulation. As in the K-theory 
simulation with surface stress removed, removal of surface stress strengthens the 
wind, geostrophic shearing deformation, and vertical motion without significantly 
effecting the strength of the cross front temperature gradient. This paradox will be 
addressed at some length below in Chapter VI. 

Stronger wind also results in stronger surface fluxes, which are fron- 

tolytic. 

In terms of both structijue and circulation this front is very similar to 
that in the previous simulation. 

6. Level 2 Simulation 

The level 2 simulation involves replacing prognostic expressions for all 
of the second-order terms with diagnostic expressions. An expression for turbulent 
kinetic energy is obtained by setting the LHS to zero and neglecting turbulent 
diffusion of turbulent kinetic energy (the first term on the RHS of equation (4.26)), 
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thereby obtaining a balance between production and dissipation. After some re¬ 
arrangement, the resulting equation has three terms on the RHS; two terms for 
shear production and one for buoyant production. 

The surface pressure and near-surface potential temperature at hours 84, 
* 96, and 108 (not shown) are very similar to those for the full physics simulation. 

A cross section of potential temperature for hour 108 is shown in Fig. 5.77. Note 
that the strength of the front and boundary layer structure in the warm sector are 
nearly identical to the full physics simulation (Fig. 5.57) but that the depth of the 
mixed layer within the frontal zone and in the cold air is much less than in the full 
physics case. Cross sections of the front-normal and along-front wind components 
(not shown) are consistent with this picture; the distributions are quite similar to the 
full physics results but features near the surface have a smaller vertical scale. The 
most important difference between these two simulations is in the vertical motion 
(Fig. 5.78). The largest ascent is less than half of that in the full physics simulation 
(5|i b s“' vs. 11 for the full physics case shown if Fig. 5.61). The vertical motion 
field is also relatively disorganized in this simulation. 

The results of this simulation show that removing turbulent mixing of 
turbulent kinetic energy results in a much shallower boundary layer and a weak, 
disorganized vertical motion field. In the full physics (level 3) simulation, turbulent 
mixing of turbulent kinetic energy resulted in the transport of turbulent kinetic en¬ 
ergy into the vicinity of the inversion, thus partially countering buoyant destruction 
due to stable stratification. This enhanced entrainment and deepening of the bound¬ 
ary layer. In the level 2 simulation, diffusion of turbulent kinetic energy is neglected 
and the boundary layer is not as deep. This reduces the depth over which turbulent 
mixing of momentum occurs which, in turn, reduces integrated mass convergence 
and vertical motion. 


73 





A level 2.5 8imulatk>n waa also performed. This involves replacing prog¬ 
nostic expressions for all of the second-order terms except turbulent kinetic energy 
with diagnostic expressions. Prognostic expressions for turbulent kinetic energy and 
the mean quantities are retained. This greatly simplifies the numerical scheme and 
level 2.5 models have gained considerable popularity. The scheme has, however, 
been criticized by Helfand and Labraga (1986) in that it may produce unphysicai 
results in situations in which the turbulence level is changing rapidly (for example, 
frontogenesis) due to neglect of the time tendency in other second-order quwti- 
ties (variances and covariances). The results of the level 2.5 simulation are nearly 
identical to those for the full physics run. This indicates that, in the present im¬ 
plementation, the neglect of the time tendency does not seriously compromise the 
results. Inclusion of moist thermodynamics may alter this conclusion. 

7. Summary 

Use of the second-order closure simulation produced a depiction of frontal 
features which was much mote realistic than the AI case and somewhat more realistic 
than the K-theory case. The numerical results are summarized in Table 5.1. Note 
that values calculated for the strong cold front from the STREX field program 
documented by Fleagle et al. (1988) are included in the last row for comparison. 
The magnitude of the potential temperature gradient across both the warm and cold 
fronts was much larger than in the K-theory simulation. The convergence, vorticity, 
and vertical motion were all much stronger than in either of the other two cases. 
The results compared favorably with the observations of maritime fronts from the 
STREX field program. 

At hour 108, the convergence and vertical motion are both three times as 
large as in the K theory simulation. The vorticity is more than four times stronger. 
The gradient in potential temperature is twice as large. Both convergence and 
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vorticity are larger than in the AI case (vorticity is more that twice as strong and 
vertical motion is seven times as large). 

The primary difference between the K-theory and second order closure 
simulations is the depth of the post-frontal mixed layer. The K-theory scheme is 
much more active in terms of turbulent mixing of heat and entrainment, both of 
which are frontolytic. This is due to differences in the length scale formulation. 
Although both schemes use the Blackadar (1962) formulation for the length scale 
(4.15), the K-theory scheme uses a constant value in the asymptotic limit (350 m for 
the eddy heat coefficient). The asymptotic limit in the second order closure scheme 
is the ratio of the first to the zeroth order moments of the profile of the turbulence 
and is thus dependent on the vertical distribution of turbulence. This results in a 
much larger mixing length above the inversion in the K theory scheme than in the 
second-order scheme. The manifestation of this difference in length scales is in the 
differing responses of the two turbulence parameterizations to the presence of an 
inversion. Examination of the vertical profiles of the eddy coefficients in the cold 
air near the front where the eddy coefficients are largest shows a relatively gradual 
decrease in the magnitude of the eddy coefficients above the inversion in the K-theory 
parameterization 1 order of magnitude (100 mb)~^) and a much sharper decrease 
with elevation in the second-order scheme (four orders of magnitude between the 
inversion and the grid point above). Previous experience with the two schemes 
is consistent with this behavior; the second-order closure scheme strengthens the 
inversion with time during entrainment more enthusiastically and realistically than 
does the K-theory scheme. Strengthening of the inversion inhibits further mass 
entrainment and deepening of the boundary layer in the second-order closure scheme. 

The physical interpretation of the differences between the second-order 
closure and K-theory schemes arises from the concept that the length scale is, in 
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some sense, related to the scale of the parameterized eddies responsible for turbulent 
mixing. In the second-order closure scheme, the scale of the parameterized eddies 
depends on the vertical distribution of turbulent kinetic energy in each vertical 
column, while in the K-theory scheme, the scale of the eddies has a constant value 
(in regions where the asymptotic limit ^plies). 

Other reasons for differences include different surface layer parameteriza¬ 
tion and stronger vortidty and convergence. Experiments were performed in order 
to determine the effect of the surface layer parameterizations on the simulations (the 
surface layer parameterizations are discussed in Chapter IV). These experiments in¬ 
dicated that the scheme used in the second-order closure simulations (the Liu et al. 
(1979) scheme) responds rapidly to changes in surface layer stability. Differences 
in stability over small spatial scales can result in large horizontal gradients in drag 
coefficients for momentum and heat. The Louis (1979) scheme used in the K-theory 
simulation does not respond as rapidly. This behavior is clear from the distributions 
of stirface heat flux discussed in Section 4 above and shown in Figs. 5.66 and 5.39. 
The larger near surface wind speed due to stronger momentum mixing in the second- 
order closure case is also partially due to use of the Liu et al. scheme. The results 
of the additional experiments performed indicate, however, that the differences in 
the surface layer parameterizations have a smaller impact on the simulations than 
do the differences in the eddy coefficients discussed above. Larger cyclonic vorticity 
in the second-order closure simulation is a consequence of stronger convergence and, 
to a lesser extent, vertical motion. Flrontogenesis in vorticity forced by vertical mo¬ 
tion is discussed in the frontal development section for the full physics second-order 
closure simulation above. Stronger convergence results from the stronger frontal 
temperature gradient and upward vertical motion. 
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VI. DISCUSSION 


A. INTRODUCTION 

The results of the adiabatic and nearly invisdd simulation demonstrate that, 
in the absence of boundary layer physical processes, a growing baroclinic wave will 
generate very strong warm and cold fronts near the surface. Relative vorticity 
and convergence within the frontal zone are large, as is the horizontal potential 
temperature gradient across the front. These simulations, however, are unrealistic 
in several ways. For example, vertical motion and the circulation about the front 
are extremely weak. 

Use of the K'theory parameterization resulted in a more realistic frontal struc¬ 
ture than in the AI case. The vertical circulation was much stronger, in spite of the 
fact that convergence was not as strong. This was due to the action of parameterized 
BL processes. 

Simulations in which various physical processes were excluded indicated that 
surface heat flux was most important in developing the frontal structure. Although 
the direct effect of surface heat flux was frontolytic, the indirect effect through 
modulation of BL structure was critical in producing realistic frontal features. Both 
convergence and vorticity in the frontal zone were smaller than those in the AI case. 
When surface heat and momentum fluxes were removed from the simulation, con¬ 
vergence and vorticity were nearly doubled and were then comparable in magnitude 
to the AI case. 

Use of the second-order closiire parameterization produced a depiction of 
frontal features which was much more realistic than the AI case and somewhat 
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more realistic than the K-theory case. The magnitude of the potential temperature 
gradient across both the warm and cold fronts was much larger than in the K-theory 
simulation. Convergence, vorticity, and vertical motion were all much stronger than 
in either of the other two cases. The results compared favorably with observations 
of maritime fronts from the STREX field program. The results of this simulation 
compared more favorably with observations than results of the K-theory simulation 
because of the more realistic treatment of turbulent mixing in the second order 
scheme due to use of a physically-based length scale formulation. 

Given the results of these simulations, several boundary layer processes im¬ 
portant in frontogenesis can be identified. The primary physical processes operative 
in simulations including boundary layer physics are surface fluxes of heat and mo¬ 
mentum and turbulent diffusion of heat and momentum. These four processes are 
manifested in many different ways and each of them gives rise to nonlinearities and 
feedbacks which hinder gaining a complete understanding of their impacts on fronto¬ 
genesis. Some of these processes are important in directly strengthening the frontal 
temperature gradient while others impact frontogenesis through modification of the 
secondary transverse circulation, which is crucial in both genesis and maintenance 
of the frontal zone. Some of these processes are, of course, important in both areas. 

B. PROCESSES DIRECTLY INFLUENCING FRONTAL 
GRADIENTS 

Turning first to surface heat flux, it is obvious that the direct, first order 
effect is on the frontal temperature gradient and is frontolytic: heat flux divergence 
associated with the upward surface flux in the cold air causes wanning, driving 
the temperature in the cold air toward that of the warm sector. Surface heat flux 
is partially responsible for weakening the front in simulations with boundary layer 
parameterizations as compared to the adiabatic and inviscid simulation. In the 
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simulation using the K-theory parameterization with surface heat flux excluded, the 
front was much stronger. Diabatic frontolysis in both the K-theory and the second- 
order closure simulations had a relative maximiun at the stirface at the leading edge 
of the front due to surface heat flux (a secondary maximum was located well behind 
the front in the cold air coincident with the maximum surface heat flux). 

In contrast to surface heat flux, the direct, first order effect of surface drag is 
on the frontal temperature gradient and is firontogenetic: frictionally-induced con¬ 
vergence at the front contracts the near surface temperature gradient. Howevo-, 
convergence not only strengthens the frontal temperature gradient directly but also 
increases cyclonic relative vorticity due to vortex stretching. The cyclonic circula¬ 
tion implied by increasing cyclonic vorticity acting on the along-front temperature 
gradient constitutes a frontogenetic (frontolytic) influence for cold (warm) fronts 
(Gidel, 1978). Thus, increased convergence has both a direct and an indirect fron¬ 
togenetic effect on cold fronts. In the adiabatic and nearly inviscid simulation, the 
bulk of the frontolysis along the warm front resulted from this mechanism. Surface 
drag also has a second-order effect which is frontolytic. As discussed by Branscome 
et al. (1989), inclusion of surface drag reduces growth in the amplitude of the 
baroclinic wave. This, in turn, reduces the strength of the geostrophic shearing de¬ 
formation. The results of the present study are consistent with those of Branscome 
et al. (1989). 

Vertical turbulent mixing of heat has several more subtle effects on frontogene- 
sis. Many of these are related to entrainment.* In the case in which the atmosphere 

certain amount of caution is necessary in identifying entrainment with turbulent mixing. 
Although entrainment results from turbulent mixing and generally increases in intensity as the 
turbulence level increases, the effects of entrainment are not entirely consistent with mixing, which 
is regarded as a diffusive process. For example, the effect of entrainment is to increase the vertical 
gradient in potential temperature at the top of the BL. However, the net effect of entrainment 
is to warm and dry the BL by combining air from the BL with warmer, drier air from the free 
atmosphere. 
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is stable aloft and at all levels in the warm sector and unstable behind the cold front, 
the growing BL behind the cold front will entrain warm air while the warm sector 
will remain quiescent. Thus, vertical mixing will be frontolytic as the cold air will 
be warmed by entrainment. This effect acts in concert with warming of the cold air 
by surface fluxes (since the entrainment flux is directly proportional to surface heat 
flux in a free convective regime). This process can also result in frontogenesis. In 
the second-order closure simulation, diabatic frontogenesis occurs within the frontal 
zone near 850 mb at hours 96 and 108. This location is near the top of the mixed 
layer where entrainment of high potential temperature air leads to warming. The 
inversion is stronger behind the front than it is in the warm sector. This distribution 
of inversion strength indicates that there is greater warming due to entrainment in 
the warm sector, which constitutes a frontogenetic effect. 

Another aspect of surface heat flux and turbulent mixing of heat is in control¬ 
ling the static stability and boundary layer depth. The importance of the control of 
boundary layer depth will be discussed below. 

An intriguing aspect of the frontal potential temperature fields shown in the 
previous chapter for the K-theory and second-order closure simulations is that the 
frontal zone exhibits no tilt in the vertical over the lowest 100-200 mb. This is in 
contrast to the adiabatic and nearly inviscid simulations (and many observational 
studies) in which the frontal zone has a pronounced vertical tilt. The well-known 
formula of Margules also mandates that the slope of a front be nonzero. The sim¬ 
ulated frontal structure is, however, consistent with the structure of marine cold 
fronts observed over the NE Pacific ocean during the Storm IVansfer and Response 
Experiment (STREX) (Fleagle et al., 1988). In order for simulated fronts to have 
a substantial tilt, the BL near the front on the cold side would have to be very 
shallow. This region, however, is very statically unstable; in both the simulation 
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and the STREX data, upward surface heat flux behind the front is nearly 200 W 
m~^. Given this strong instability, a deep mixed layer would be expected. Thus, 
it appears to be unreasonable to expect maritime cold fronts to exhibit the vertical 
tilt commonly associated with continental cold fronts. 

One effect of the the ageostrophic vertical circulation is to cause the front 
to tilt (Bluestein, 1986). In the following discussion, we will refer to this as the 
‘‘sloping effect” in order to avoid confusion with the vertical deformation term of 
the frontogenetical function. In the upper portion of the BL, warm air is advected 
toward the cold air by the horizontal branch of the circulation while, in the lower 
BL, cold air is advected toward the warm air, resulting in a tilted front. From the 
Sawyer-Elliassen equation, we can infer (following Bluestein) that, if inertial stability 
(absolute vorticity in the term involving the vertical derivative of the ageostrophic 
momentum) is large compared to static stability (in the term involving the horizontal 
derivative of the vertical motion), vertical motion dominates and horizontal motion is 
inhibited; if absolute vorticity is small compared to static stability, horizontal motion 
dominates and vertical motion is inhibited. In the second-order closure simulations, 
the static stability is near zero while the relative vorticity is approximately 0.6 x 10~^ 
s~^ so that horizontal motion is inhibited and vertical motion dominates. Thus, the 
sloping effect of the vertical circulation will be reduced. This provides a dynamic 
mechanism for production of vertical fronts in the presence of strong instability. 
In a study of a continental cold front with a substantial vertical tilt, Ogura and 
Portis (1982) found the maximiim relative vorticity to be 1 x 10“^ s”* . The static 
stability (the Brunt-Vaisala frequency) in this case can be crudely estimated from a 
cross section shown in their Fig. 17 as ^ 1.5 x 10”^ s~*. Therefore, in this continental 
cold front exhibiting substantial vertical tilt, the static stability is large compared 
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to the relative vorticity and horizontal motion dominates. Thus, the sloping effect 
of the vertical circulation is enhanced. 

It is ircmic that the primary criticism of quasi-geostrophic fronts is that they 
have no tilt in the vertical (Stone, 1966) pven the present results and those of 
recent observational studies showing that maritime fronts ^parently are vertical 
at low levds. There are, however, other characteristics of quasi-geostrophic fronts, 
which are referred to as ‘‘pseudo-fronts” by Williams (1967), not shared by fronts 
depicted in the current simulations. For example, the vorticity distribution of fronts 
in the present study is more realistic than b the case with quasi-geostrophic pseudo¬ 
fronts which have regions of strong negative relative vorticity adjacent to the front 
on the cold side. In addition, pseudo-fronts exhibit extremes of static stability 
with large positive vertical gradients in potential temperature in the cold lur and 
negative gradients in the warm sector. These unrealbtic features greatly enhance 
the frontolytic effect of adiabatic temperature changes due to vertical motion. 

C. PROCESSES INFLUENCING THE TRANSVERSE VERTICAL 
CIRCULATION 

The importance of the thermally direct ageostrophic transverse secondary cir¬ 
culation in both frontogenesis and in the maintenance of mature frontal zones can¬ 
not be over-emphasized. Convergence at the surface under the rising branch of 
the circulation acts frontogenetically, and augments the effects of geostrophic de¬ 
formation while vertical deformation increases the cross-front temperature gradient. 
Thus, synoptic-scale, geostrophic deformation b regarded as initiating frontogene¬ 
sis, whereas the completion and maintenance of the front are due to the vertical 
circulation (Eliassen, 1959). 

The secondary ageostrophic circulation b driven by the dynamics of frontogen- 
esb. It acts to restore thermal wind balance by 1) adiabatic temperature changes due 





to vertical motioii (which reduce the crose>fr(»t teaq>ereture gradioit) (Blueetem, 
1986) and 2) increasing the vutical shear much as the setxmdaiy vertical circulation 
associated with a growing baroclinic wave a^ to keq> vortidty changes geostrophic 
and temperature changes hydrostatic. The geostrophic vertical shear is enhanced 
due to the frontal baroclinic zone (thermal wind). Vertical shear in the along-front 
wind increases towards the geostrophic value by Coriolis turning of the front-normal 
ageostrophic wind: 

^ = -/(kxV.). 

The ageostrophic wind is a response to the increased slope of pressure surfaces. 

Inclusion of BL physics augments the strength of the ageostrophic circula¬ 
tion and thereby contributes to frontogenesis. Surface drag enhances low-level con¬ 
vergence. Convergence acts frontogenetically in that it enhances the strength of 
the ascending branch of the secondary ageostrophic circulation. Vertical deforma¬ 
tion associated with the drculation increases near-surface baroclinicity and vorticity 
within the frontal zone. Note also that friction opposes its own dissipative effect on 
vorticity.* 

In both the K-theory and second-order dosure simulations, vertical deforma¬ 
tion was the primary source of frontogenesis while, in the adiabatic and nearly 
invisdd simulation, the contribution from tilting was much smaller. This resulted 
from small upward motion in the adiabatic and inviscid case. As discussed d>ove 
in Chapter V, while the pattern of vertical motion was similar in all three of these 
simulations, the magnitude was much smaller in the adiabatic and inviscid case. 

*Wuids an convergent (divergent) in ngions of cyclonic (anti-cyclonic) relative vorticity due to 
nirface friction. Convergmoe enhances nUtive vorticity. Thus, friction opposes its own dissipative 
effect on vorticity throng the divergence term (Carlson, 1991). This idea was alluded to by KA82 
when they stated that "cmivergence in the frictionally driven PBL planetary boundary lay«] 
inflow... counteracts the dissipation d vorticity by PBL friction...”. 
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An important result from the K-theory amulation is that, while the conver¬ 
gence is only about half the value in the A1 case, the updraft is more than twice as 
strong. This results from convergence extending through a deep layer. Thus, the 
turbulent mixing in concert with the surface fluxes rurganises convergence in such a 
way as to produce a mudi more compact and vigorous vertical drculatkm. This is 
largely accomplished by mixing momentum near the surface so that the horizontal 
shear exists through a relatively deep layn. In the second-mder closure simulation, 
both convergence and vertical motion are stronger than in the K-theory case. When 
the effect of surface drag is removed, the vertical motion is much stronger. This 
leads to the paradoxical conclusion that the transverse ageostrophic circulation is 
stronger in the absence of surface drag. Closer examination of the results, however, 
indicates that convergence at and near the surface is nearly the same as in the full 
physics simulation. The effect of removing surface drag is evident in the wind field 
(particularly in the front-normal componoit). The wind is stronger but the strong 
wind extends almost to the surface. Thus, although the wind is much stronger, the 
convergence is not. This is an encouraging result in that it verifies the importance 
of surface drag in enhancing low level convergence. The explanation for enhanced 
vertical motion lies in mixing of momentum in this case also. Strong convergence 
extends all the way to the surface as well as existing through a relatively deep layer, 
giving rise to a large value of integrated mass convergence. In the full physics sim¬ 
ulation, wind speed and convergence decrease towards the surface due to surface 
drag. These results are summarized in Table 6.1. 

The results of the level 2 simulation are consistent with the argument that tur¬ 
bulent mixing of momentum b critical in strengthening the transverse ageostrophic 
circulation as well. When turbulent mixing of turbulent kinetic energy is removed, 
the vertical motion b less than half as strong as in the full physics simulation. 
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In ■iiininaiy, theie nrete indicatioiu tm to the inqMrUnoe ot turbulent mi x i n g of 
momentum in the K-theory simulation. The nxne realistic treatment of turbulence 
in the second order closure simulation clarified this issue. The results are consistent 
with the idea that frictional depletion oi momentum at the surface combined with 
turbulent mixing will cause a momentum deficit throughout the BL. This in itself 
creates ageostrophic flow since the Coriolis force (fv) will be too small rdative 
to {l/p)dp/dz (see Keyser and Anthes, 1982). The important point is that this 
ageostrophic flow will be directed towards the front, enhancing convergence over a 
deep layer. The integrated mass convergence then results in strong upward motion. 
This also shows that, while the direct effect of surface heat flux and turbulent mixing 
of heat is frontolytic, the role that heating plays in modulating boundary layer depth 
is important in frontogenesis. 

These results are, to a certain extent, anticipated by the work of Nuss (1989). 
In simulations of explosive marine cydogenesis, he found that the prefrontal updraft 
was driven by frictionally induced boundary layer convergence and the distribution 
and intensity of this convergence were influenced by stability. F^hermore, the 
surface wind stress, which is a measure of the downward momentum fluxes in the 
BL, is also influenced by boundary layer stratification changes. The surface wind 
stress was found to be largest in unstable regions. Thus, we may conclude based on 
the results of Nuss that the prefrontal updraft will be enhanced in unstable r^ons 
due to increased downward momentum fluxes — a conclusion consistent with that 
of the present study and with the results of both the K-theory and second-order 
closure simulations. 

This discussion has a bearing on the arguments presented by Koch (1984), 
Pinkerton (1978), Dorian et al. (1989), and Reeder (1986). These authors discuss 
scale contraction and intensification of the vertical circulation in the presence of 
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horiKmtolly differential heat flux (with the strongest heat flux in the warm sector). 
Reeder explains this based on higher order terms in the divergence equation leading 
to enhanced convergence in areas experiendng convergence in an insulated surface 
simulation. This argument has not gained widespread acceptance. In Reeder’s 
* simulation, a cold front moved onshore and encountered a deep mixed layer in the 
warm sector caused by strong upward heat flux from ground heated by insolation. 
In other studies, the warm sector was he^ed by horizontally differential surface heat 
flux arising from other causes. In the present study, there is also an intensification of 
the vertical circulation in the presence of surface heat flux. The heat flux, however, 
is largest in the cold air behind the front and constitutes a frontolytic effect. This 
argument is supported by comparison of results of the K-theory to the second-order 
closure simulation. Recall that, in the K-theory simulation, the horizontal gradient 
in surface heat flux was extremely weak while, in the second order closure simulation, 
the difference across the front was in excess of 120 W m~’. The upward vertical 
motion in the second-order closure simulation was three times as large as in the 
K-theory simulation. A common theme between the present results and published 
work is the relatively deep mixed layer in the warm sector. Thus, in all of these 
cases, vertical turbulent mixing of momentum was taking place in and ahead of the 
frontal zone. In light of the present results, it is likely that mixing of momentum in 
the deep mixed layer in and ahead of the frontal zone may explain intensification of 
the vertical transverse drculation in the presence of horizontally differential surface 
heat flux. 

The adiabatic and invisdd, K-theory, and second-order closure parameteriza- 
tions comprise a hierarchy in levels of sophistication in treatment of vertical tur¬ 
bulent mixing. The strength of the vertical circulation increases with increasing 
sophistication. In the adiabatic and invisdd simulation, there is small vertical mo- 
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tkm while, in the leoondorder ckMuie simttlntion, the accent exceeds that in come 
published (^Mervatkmal studies. 

Keyser and Anthes (1982) state that ^...increasing the vertical resolution in 
the PBL results in improved vertical structure of the... u jet... Clearly, increasing 
the vertical resolution is not adequate. The adiabatic and inviscid, K-theory, and 
second-order closure simulations were all run with the same vertical resolution but 
only in the second-order closure case is there strong vertical motion resulting from 
more realistic treatment of turbulent mixing. 
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VIL CONCLUSIONS 


The primary conclusioiu of the present study are as follows: 

1. A hydrostatic primitive equation three-dimensional model simulated the growth 
of a baroclinic wave from a small amplitude perturbation on a barodinically 
unstable flow to a deep low with strong circulation and warm and cold fronts. 
Lower boundary conditions were specified in a manner consistent with oceanic 
cyclogenesis/frontogenesis. Analyses of the structure and dynamics of the 
growing baroclinic wave at 72 h showed its behavior to be consistent with 
observations and traditional conceptual models. It may thus be concluded 
that the approach used in this investigation (the use of a numerical model to 
simulate the growth of a baroclinic wave and development of the associated 
warm and cold fronts) is sound. 

2. The adiabatic and nearly invisdd simulation produced intense fronts with 
strong convergence and vortidty but with an unrealistically weak vertical 
transverse circulation. The primary process leading to frontogenesis in this 
case was geostrophic shearing deformation. 

3. The K-theory simulation produced a somewhat more realistic depiction of 
frontal structure, although the vortidty and convergence were weaker than 
in the adiabatic and invisdd case. Several boundary layer physical processes 
played a role in producing more realistic results. Specifically, surface heat 
and momentum fluxes combined with turbulent mixing of heat promoted the 
development of a deep, well-mixed boundary layer. Surface drag also reduced 




near-surface wind speed and produced frictional convergence. IVirbulent mix¬ 
ing of momentum combined with frictional convergence lead to the devebp- 
ment of a stronger transverse secondary circulation than in the AI case. The 
strength of the front was greatly reduced by surface heat flux. 

4. Several numerical experiments were performed in which various physical pro¬ 
cesses were excluded from the K-theory simulation. The results suggested that 
surface heat flux was the most significant process in producing the results of 
the full physics simulation. 

5. The second-order closure simulation produced stronger vorticity, convergence, 
and vertical motion than either the K-theory or the adiabatic and invisdd 
simulation. The front was twice as strong as in the K-theory case but not 
nearly as strong as in the adiabatic mid invisdd case. Results of this simula¬ 
tion compared favorably with observations of maritime fronts from the Storm 
Transfer and Response Experiment (STREX). 

6. Results of numerical experiments in which various physical processes were 
excluded indicated that both surface heat and momentum fluxes were required 
to reproduce the results of the full physics simulation. 

7. Examination of the results indicated that the use of a length scale based on 
the vertical distribution of turbulent kinetic energy in the second-order closure 
scheme resulted in a more realistic parameterization of turbulence than the K- 
theory scheme, in which the asymptotic limit of the length scale b a constant. 
This length scale prescription resulted in a sharp decrease in eddy coeffidents 
with elevation above the inversion in the second order closure case and more 
gradual decrease in the K-theory case. This K-theory formulation resulted 
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in enhanced entrainment and a mixed layer whi^ waa too deq> and a front 
which was too weak as compared to observations from STREX. 

8. In both the K-theory and second order closure simulations, the fronts did 
not tilt in the vertical over the lowest 100-200 mb. This structure is broadly 
consistent with observations of fronts over the NE Pacific ocean and apparently 
results from strong upward surface heat flux in the cold air behind the front. 
Dynamical arguments involving the transverse vertical circulation suggest that 
vertical fronts may be expected in rejpons of neutral static stability. 

9. Analyses of frontogenetic forcing show that vertical deformation was most im¬ 
portant in both the K-theory and second-order closure simulations. Stretching 
and shearing deformation played small but significant roles in both simulations. 
In the adiabatic and invisdd simulation, stretching and shearing deformation 
were important along different portions of the front. Vertical deformation was 
very small in the vicinity of the fronts. 

In the K-theoiy and second-order closure simulations, frontolysis along both 
the warm and cold fronts was due to diabatic frontogenesis, which had a rel¬ 
ative maximum at the surface at the leading edge of the front. In aU three 
simulations, frontolysis along the warm front was due to shearing deforma¬ 
tion, which quantifies the impact of cyclonic circulations on the along-front 
temperature gradient. 

10. In both the K-theory and second-order closure simulations, vertical turbulent 
mixing of momentum enhanced vertical motion and, therefore, frontogenesb. 
In the K-theory simulation, upward motion ahead of the frront in the warm 
sector was stronger than in the adiabatic and invbdd case in spite of the fact 
that the convergence was not as large. In the second-order closure simulation. 
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cmivergence and vertical motion were both stronger than in the K-theory sim¬ 
ulation. In the second-order dosure simulation with surface drag removed, 
the convergence near the surface was only slightly stronger than in the full 
physics simulation but the vertical motion was much stronger due to vertical 
» turbulent mixing of momentum. 

11. The results suggest the following scenario for midlatitude maritime frontoge- 
nesis: 

(a) Geostrophic shearing deformation assodated with a growing baroclinic 
wave begins to develop a strong baroclinic zone on the northeast side 
of the low with a weaker zone extending south and west at 72 h. At 
this time, there is large upward motion due to strong convergence in the 
vidnity of the strong baroclinic zone with little descent and no significant 
vertical motion assodated with the weaker zone. Boundary layer physical 
processes do not play an important role at this time. 

(b) Twelve hours later, the northern baroclinic zone, now recognizable as 
a warm front, extends over the low from west to east and southeast in 
advance of the low. The second baroclinic zone, now recognizable as a 
cold front, extends southwestward from the south side of the low. The 
two fronts do not intersect at the low but are separated at the frontal 
fracture. The thermal gradients across both zones are nearly the same as 
at 72 h. The upward motion in the vicinity of the warm front is also the 
same as at 72 h. There is weak upward motion and no descent assodated 
with the southern branch of the warm front. There is now a strong, 
thermally direct drculation about the cold front and tilting frontogenesis 
begins to accelerate contraction of the thermal gradient at this location. 
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Weak tilting firontogenesis along the warm front ia cancded by shearing 
deformation. 

(c) At 96 h, the cold front has strengthened dramatically while the northern 
and southern portions of the warm front have weakened slightly. The 
direct circulation about the cold front b now stronger than at any other 
time. The boundary layer in thb repon is approximately 150 mb deep 
and vigorous mixing of momentum is occurring along with frictional con¬ 
vergence, augmenting the strength of the updraft ahead of the front. 
Strong upward surface heat flux and entrainment are occurring. Both 
of these processes are frontolytic. Although it has weakened over the 
last 12 hours, there is still significant upward motion associated with the 
northern branch of the warm front and little vertical motion along the 
southern branch. Tilting frontogenesb is still quite strong along the cold 
front and small along the southem portion of the warm front. The con¬ 
tribution due to tilting along the warm front b canceled by frontolysb in 
both shearing and stretching deformation. 

The warm and cold fronts now intersect and an occlusion is developing 
to the north of this intersection. 

(d) At 108 h, the strength of the cold front is unchanged. Tilting fronto- 
genesis, which is still strong, b balanced by diabatic frontolysis with a 
smaU contribution from the shearing deformation. Boundary layer depth 
is unchanged. The northern and southern portions of the warm front are 
again slightly weaker than at 96 h. FVontolysb in the shearing deforma¬ 
tion now dominates the total adiabatic frontogenesb along the southern 
portion of the warm front. 
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The occlusion is now a prominent feature of the simulation. Horismital 
cross sections spanning the domain at different latitudes clearly show 
warm air being forced aloft at higher latitudes. 

(e) At 120 h, all of the fronts are weaker than at 108 h. The occlusion is now 
the most prominent feature. 

In the near future, this work should be extended in simulations using higher 
horizontal resolution. This should be possible with currently available computers.^ 
Some work in this direction has been performed unng nested grids (Gall et al., 1987). 
In time, it may be possible to understand more fully the specific roles fulfilled by 
turbulent mixing in frontogenesis through the use of large eddy simulation (LE)S) 
of fronts and baroclinic waves. However, this work must await the development of 
more powerful computers; current LES implementations require hundreds of hours 
of cpu time for domains on the order of lO's of km-simulation of fronts and baroclinic 
waves requires domain sizes two orders of magnitude larger. 

Another area in which further work is required is in documenting the char¬ 
acteristics of the westward extension of the warm front (WEWF) in the Shapiro 
and Keyser (1990) conceptual model. Some effort in this direction was made in the 
present study (see Appendix). The results, although incomplete, are encouraging 
but not entirely consistent with that model. For example, Shapiro and K^ser do 
not include an occlusion in their conceptual model This was not the case in either 
the adiabatic and invisdd simulation or the second-order closure simulation. Many 
of the results concerning the seclusion and the WEWF are consistent with more 
recent publications by Shapiro and his co-workers (e.g., Neiman and Shapiro, 1993; 
Neiman et al., 1993; Kuo et al., 1992). FWther study in also required on the impact 

^Tbe present study was performed on a CY205 computer; large amounts of computer time and 
the entire central memory were requited for the seconc^rder closure simulations. 
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of differing sea surface temperature distributions on the WEWF and seclusion and 
interaction of the WEWF with the occlusion. 

In the present study, two different surface layer parameterization sdiemes were 
used. The more sophisticated scheme used in the second-order closure simulations 
produced larger horizontal and vertical gradients in fluxes of heat and momentum. 
In the future, additional surface layer elements should be investigated. In particular, 
the effects of wave age and sea spray should be included. In the vicinity of a maritime 
front, unsteadiness in the wind field due to changes in fetch and wind speed and 
direction could significantly impact the evolution of roughness elements. In addition, 
sea spray could have a large effect on the drag coefficient for heat and moisture. 


95 




[THIS PAGE INTENTIONALLY LEFT BLANK] 


96 








lkbl«4.1: MODEL CHARACTERISTICS 



Value 

Donudn Size (x) 

2960 km 

Domam Size (y) 

4000 km 

Ax 

40 km 

Ay 

40 km 

Ax 

6 

At 

120 seconds 

/ 

1.031 X 10“** 

Jet Max 

40 m/s at 200 mb 

Meridional Jet Profile 


Perturbation Amplitude 

1 m/s 

No, of Gravity Modes 

1 m/s 

Horiz. Diffusion Coeff. 

-10* 





Tkble4.a: VERTICAL STRUCTURE 
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Tkble 6.1: FULL PHYSICS VS. NO DRAG 




No Drag 

Percentage 

Change 

u» in h/s) 

11.11 

13.89 

20% 

Normal Component (m/s) 


14.00 

55% 

Divergence (lO'^/s) 

1.10 

1.20 
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Figure 5.2: Meridional cross sections of the initial a) potential temperature (K) 
and b) zonal wind component (m s"*) from the surface to 10 mb and from the 
northern boundary to the southern boundary. 
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Figure 5.2: Continued. 
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Figure 5.3; Time series of minimum surface pressure (mb) for adiabatic and in- 
visdd simulation (AI), K-theory simulation (KTH), and second-order closure simu¬ 
lation (20C). 
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Figure 5.4: a) surface pressure (mb) and b) near-surface potential temperature 
(K) fields at hour 72. 
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Figure 5.4: Continued. 
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Figure 5,5: Continued. 
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Figure 5.5: Continued. 
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Figure 5.6: Near-surface potential temperature (K) fields at hour 84 (AI). 
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5.8: Near-surface potential temperature (K) fields at hour 108 (AI) 
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Figure 5.10: Continued. 
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Figure 5.11: Vertical cross section of potential temperature (K) at hour 72 FV(»n 
the surface to 500 mb and firom 0 to 1320 km. The plane of the cross section, 
denoted by AB, is shown in Fig. 4b (AI). 
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Figure 5.12: Cross section of potmitial temperature (K) at hour 84. The plane 
of the cross section, denoted by CD, is shown in Fig. 6 (AI). 
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Figure 5.13: Cross section of potential temperature (K) at hour 96. The plane 
of the cross section, denoted by EF, is shown in Fig. 7 (AI). 


118 













Q 


Oistanctt (km) 


40 km 


H 


Figure 5.15: Cross sections of a) along-front wind component and b) front-normal 
wind component (m s~^) at hour 108 (AI). 
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Figure 5.16: Cross section of divergence (xlO~^ s“*) at hour 108 (AI). 
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Figure 5.17: Cross section of vertical motion (/i b s~^) at hour 108 (AI). 
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Figure 5.18: Cross section of relative vorticity (xlO~^ s~*) at hour 108. The 
plane of the cross section, denoted by GH, is shown in Fig. 9 (AI). 
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Figure 5.19: Continued. 
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Figure 5.19: Continued. 





















Figure 5.20: Total adiabatic frontogenetic forcing at hour 108 at 850 mb (AI). 
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Figure 5.21: a) cross section of total adiabatic frontogenetic forcing (K s 
at hour 108 and b) stretching deformation (AI). 
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5.23: Near-surface potential temperature (K) fields at hour 84 
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Figure 5.24: Near-surface potential temperature (K) fields at hour 96 (KTH). 
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Figure 5.25: Near-surface potential temperature (K) fields at hour 108 (KTH). 
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Figure 5.26: 850 mb relative vorticity (xl0“* s“*) at hour 108 (KTH). 


136 













0 Distance (km) x 


3000 


Figure 5.27: a) surface pressure (mb) and b) near-surface potential temperature 
(K) fields at hour 120 (KTH). 
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Figure 5.27; Continued. 













Figure 5.28: Cross section of potential temperature (K) at hour 84. The plane 
of the cross section in shown in Fig. 23 (KTH). 
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Figure 5.20: Cross section of potential temperature (K) at hour 96. The plane 
of the cross section is shown in Fig. 24 (KTH). 
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Figure 5.30: Cross section of potential temperature (K) at hour 108. The plane 
of the cross section is shown in Fig. 25 (KTH). 
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Figure 5.31: a) cross section of along-front wind component and b) front-normal 
wind component (m s~*) at hour 108 (KTH). 
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Figure 5.32: Cross section of divergence (s~^) at hour 108 (KTH). 
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Figure 5.33: Cross section of vertical motion b s~^) at hour 108 (KTH). 
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Figure 5.34: Cross section of vorticity at (x 10~^ s~^) hour 108. The plane of the 
cross section is shown in Fig. 26 (KTH). 
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Figure 5.37: Cross section of total adiabatic frontogenetic forcing at hour 96 
(KTH). 
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Figure 5.38: a) cross section of total adiabatic frontogenetic forcing at hour 108, 
b) tilting term, and c) diabatic term (KTH). 
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Figure 5.38: Continued. 


151 






















Pressure (mb) 
















Figure 5.39: Surface sensible heat flux (W at hour 108 (KTH). 
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Figure 5.41: Near*surface potential temperature (K) fields at hour 96 (KTH no 
flux). 
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Figure 5.43: a) cross section of potential temperature (K) at hour 84, b) at hour 
96, and c) at hour 108. The planes of the cross sections are shown in Figs. 40, 41, 
and 42 (KTH no flux). 


157 


Pressure (mb) 





Figure 5.43: Continued. 
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5.44: a) normalized eddy coefficient for momentum {Km) and b) normal¬ 
ly coefficient for heat (Kh) (KTH). 
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Figure 5.46: Cross section of potential temperature (K) at hour at hour 108. The 
plane of the cross section is shown in Fig. 45 (KTH no stab). 
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Figure 5.47: Time series of minimum surface pressure (mb) (20C). 
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Figure 5.48: Near-surface potential temperature (K) fields at hour 84 (20C). 
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Figure 5.49: Near-surface potential temperature (K) fields at hour 96 (20C). 
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Figure 5.50: Near-surface potential temperature (K) fields at hour 108 (20C). 
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Figure 5.51: 850 mb relative vortidty (xlO"® s“') at hour 108 (20C). 
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Figure 5.52: a) surface pressure (mb) and b) near-surface potential temperature 
(K) fields at 120 hour (20C). 
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Figure 5.52: Continued. 
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Figure 5.53: Cross section of potential temperature (K) at hour 84. The plane 
of the cross section is shown in Fig. 48 (20C). 
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Figure 5.54: Cross section of potential temperature (K) at hour 96. The plane 
of the cross section in shown in Fig. 49 (20C). 
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Figure 5.55: Cross section of vortidty (xlO~^ s~*) at hour 96 (20C). 
























Figure 5.57: C^oss section of potential temperature (K) at hour 108. The plane 
of the cross section is shown in Fig. 50 (20C). 












Figure 5.58: Continued. 
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Figure 5.60: 


Cross section of divergence (xlO~^ s"^) at hour 108 (20C). 
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Figure 5.64: Cross section of total adiabatic frontogenetic forcing at hour % 
(xlO-« K m-* s->) (20C). 
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Figure 5.65: a) cron section of total adiabatic frontogenetic forcing K 

s”*) and b) tilting frontogenesis at hour 108. The plane of the cross section is 
shown in Fig. 62 (20C). 
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Figure 5.65: Continued. 
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Figure 5.69: Near-surface potential temperature (K) field at hour 108 (20C no 
flux). 
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Figure 5.70: a) cross section of potential temperature (K) at hour 84 (the plane 
of the cross section is shown in Fig. 5.48), b) at hour 96 (see Fig. 5.49), and c) at 
hour 108 (see Fig. 5.50) (20C no flux). 









Figure 5.70: Continued. 
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Figure 5.71: Near-surface potential temperature (K) field at hour 84 (20C no 
drag). 
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Figure 5.75; Cross section of divergence (xlO“^ s”*) at hour 108 (20C no drag). 
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Figure 5.76: Cross section of vertical motion (/i b s'**) at hour 108 (20C 
drag). 
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Figure 5.77: Cross section of potential temperature at hour 108 (20CL2). 
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APPENDIX 

T-BONE STRUCTURE 

In a recent paper, Shapiro and Keyser (1990) (hereinafter SK) present a ‘‘conceptual 
visualization of cydone-firontal evolution” which is ronarlcable in that it r^resents 
one of the few substantiative alteration of the Norw^ian frontal cycbne model in 
70 years. This model was originally presented by Bjerknes (1919) and Bjerknes 
and Solberg (1921 1922). The model, which enjoyed initial wide acceptance, was 
further supported by observational studies (Bjerknes and Palmen, 1937; Sanders, 
1955). In the mid 1960’s, the first weather satellites presented eloquent testimony to 
the veracity of the Norwegian model. According to Reed (1990), numerous examples 
can be found (in early satdlite imagery) in which the cloud patterns conform at each 
stage to the classical (Norwegian) scheme. Nevertheless, as Shapiro and Ke 3 rser note, 
recent observational and numerical studies of cyclogenesis and frontogenesis display 
certain features not accounted for in this model. 

The fundamental departure from the Norwegian model occurs in the second 
stage of frontal evolution. In stage 2 of the Norwepan model, the temperature wave 
has attained sufficient amplitude that the warm and cold fronts are discernible and 
the cloud shield begins to develop (Fig. A.1). In stage 2 of the SK model, as the 
temperature wave attains similar amplitude, ”frontal fracture” occurs in which the 
warm and cold fronts separate and the warm front extends to the west into the 
northerly flow west of the cyclone center (Fig. A.2). In subsequent stages in the 
Norwegian model, the fronts remain connected and an occlusion forms in the final 
stage. In the SK model, the warm front continues to wrap around the cyclone and 
a warm seclusion forms in the final stage. The absence of the occluded stage in the 
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SK model is a major dei>arture from the Norwegian model. Of equal importance, 
however, is the presence of the westward extension of the warm front (WEWF), 
which is not a feature of the Norw^ian model at all. In many cases, the baroclinic 
zone associated with the WEWF is is even stronger than that associated with the 
cold front. Both the absence of the occlusion and the presence of the WEWF 
are directly attributable to the frontal fracture occurring in stage 2. The authors 
emphasize that the essence of the frontal evolution in their model was realistically 
simulated using adiabatic numerical models. While boundary layer processes modify 
the time scales and intensity of baroclinic waves, the basic structural evolution is 
driven by the adiabatic components. 

The authors offer no explanation for the strength of the WEWF or for the 
frontal fracture, which is manifested as a loss of baroclinidty in the cold front near 
the cyclone center. The strength of the WEWF may be due to the shearing defor* 
mation; while cyclonic shear is frontolytic for warm fronts, it is frontogenetic for 
cold fronts and for the WEWF due to the sense of the along front temperature 
gradient. The loss in baroclinicity near the cyclone center can be explained using 
the results of Keyser et al. (1988). In this analytical study, an E-W band of isen- 
tropes is superimposed on an axisymmetric vortex. As this pattern evolves, strong 
frontogenesis occurs in the NE and SW quadrants, however, while the isentropes 
near the center of the vortex undergo rotation, the potential temperature gradient 
remains unaltered. This can be explained in terms of the angles between the axes of 
dilatation and the potential temperature gradient. The authors include an analysis 
of the vector frontogenetic forcing which reveab that the maximum magnitude is in 
the NE and SW quadrants while the rotational component has a relative maximum 
near the crater of circulation. 
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The results of the adiabatic and inviscid simulation in the present study axe 
not entirely consistent with the SK model. While the frontal fracture, seclusion, 
and WEWF are present, there are also strong indications that an occlusion has 
formed. Frontal fracture occurs at hour 72. The results of the adiabatic and inviscid 
simulation show that the WEWF weakens from hour 84 (Fig. A.3) to hour 120 (Fig. 
A.4). The WEWF undergoes strong frontogenesis from hour 72 to hour 84. At this 
early stage, strong convergence leads to frontogenesis. It u not clear presently why 
the WEWF weakens after hour 84. FVirther investigation will be required to clarify 
this point. The warm frontal seclusion also is displayed in Fig. A.4. This figure is 
similar to Fig. 10.19 of SK. An alternative depiction of the seclusion and WEWF is 
given in Fig. A.5, which shows the near-surface potential temperature distribution at 
120 h. Note the strength of the northern section of the WEWF, which b comparable 
to that of the cold front at this time. 

Evidence for the existence of an ocdusion can be seen in cross sections of po¬ 
tential temperature from three planes at different latitudes for hour 108 (not shown). 
The warm air is seen to rise to the north from one plane to another. The warm air 
can be delineated by the location of the 288 K isentrope. In the southernmost cross 
section, the 288 K isentrope intersects the surface at two locations 400 km apart. 
In the intermediate cross section 400 km to the north, the intersections are 200 km 
apart and, in the northernmost cross section 400 km further north, there is only one 
intersection. 

It is curious that the simulation produces both a seclusion and an occlusion. 
This suggests that the absence of an occlusion in the SK model constitutes a flaw. 
Close examination of the results suggests that the seclusion forms in the early stages 
of cyclogenesis between hour 72 and hour 96. During this time, warm air flows 
northward into the center of the low. The near surface potential temperature is 
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above 290 K at hour 84 in the center the low. Near hour 96, the frontal fracture 
I 4 >pear 8 to "heal” and the warm and cold fronts re-unite in the southeast quadrant 
of the low. At this point, the flow of warm air into the center aS the low is halted as 
the warm air begins to move aloft in the manner of the classical occlusion process. 

* The results indicate that the potential temperature in the center of the low never 
exceeds 290 K and the area enclosed by the 290 K isentrope grows smaller with time 
as the occlusion grows larger. 

It thus appears that the fourth and finAl stage of the SK model is actually an 
intermediate stage. The simulation strongly suggests that, in the final stage, the 
warm and cold fronts re-unite, sealing the seclusion and forming an occlusion. 

Published work by Shapiro and his co-workers subsequent to SK (e.g., Neiman 
and Shapiro, 1992; Keyser, 1993, personal communication) have discussed the pres¬ 
ence of occlusions at various stages of cyclogenesis. This suggests that the model as 
presented by SK has perhaps been amended to some extent to include the possibility 
of an occlusion in later stages. 

The results of the second-order closure simulation are, not surprisingly, quite 
different from those for the adiabatic and inviscid simulation. The inclusion of 
boundary layer physics has a profound effect on the simulation. Figure A.6 shows a 
cross section of potential temperature at hour 84. While the warm front, cold front, 
and WEWF are clearly present, they are much weaker in the present simulation. 
There is no indication of a seclusion near the surface. Shown in Fig. A.7 is a north- 
south cross-section of potential temperature at hour 108, corresponding to Fig. A.4 
(although there is a 12 h time difference). 

In this case, the seclusion is not as well-defined as in the AI case and does 
not extend to the surface. The southern side of the outward sbping baroclinic ring 
is not well-defined at hour 108, although there is a weak baroclinic zone extending 
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from approximately 850 mb to 600 mb delineated by the 291 and 294 K iaentropes. 
At hour 84, the seclusion is most apparent at 850 mb. FVom hour 84 to hour 108, 
the seclusion extends upward and is most well-defined at 700 mb by hour 108 (see 
Fig. A.8). This is entirely consistent with the results of Neiman and Shapiro (1993). 
In their investigation of a cyclogenesis event during ERICA, Neiman and Shi^iro 
found the first evidence of a seclusion at 850 mb. The WEWF and seclusion then 
developed upward, reaching 500 mb after 18-24 h. The relatively small scale of the 
seclusion is also consistent with Neiman and Shapiro. 

Strong evidence for the existence of an occlusion in this case is provided by 
e xamina tion of cross sections at different latitudes (not shown). The warm air can 
be delineated by the position of the 296 K isentrope. In southernmost cross section, 
the 296 K isentrope intersects the surface at two locations approximately 280 km 
apart. At an intermediate location 400 km to the north, the intersections are 50 km 
apart and, in the northern cross section 400 km further north, the 296 K isentrope 
is nowhere lower than 750 mb. 

SK present a modification to the classical Norwegian polar frontal cyclone 
which incorporates many of the features depicted in recent observational and nu¬ 
merical investigations of the life cycles of baroclinic waves and fronts. The funda¬ 
mental differences between the SK model and the Norwegian model are the absence 
of an occlusion, the presence of a seclusion, and the westward extension of the warm 
front. Both the absence of an occlusion and the WEWF are attributable to the 
frontal fracture which occurs in the early phase of the amplification of the tempera¬ 
ture wave. The authors claim that adiabatic numerical models capture the essence 
of frontal evolution in their model. 

The results of the adiabatic and inviscid simulation are, to a large extent, con¬ 
sistent with the model proposed by SK. The frontal fracture, WEWF, and seclusion 
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«re all pment in tbe linmUtion. There k abo, however, an ooclumm in the remits. 
The seclusion appears to form in the early sta^ oi cyclogenesis, prior to the for¬ 
mation of the occlusion. As the occlusion begins to form, the flow of warm air into 
the center of the low is halted. 

• The results of the second-order closure simulation are quite different. The 

warm and cold fronts and the WEWF are all much weaker in this simulation. There 
is no indication of a seclusion at the surface. A north-south cross section through 
the WEWF, however, shows some indication of a weak devated seclusion. This 
simulation also produces a well-defined occlusion. 

In conclusion, the results suggest that introduction of concepts of the frontal 
fracture and WEWF in the SK model constitute an important contribution to the 
understanding of the life cycle of baroclinic waves and fronts. The shortcoming of 
the model is that it fails to account for the last stage of the life cycle in which the 
model shows the warm and cold fronts reuniting, sealing the seclusion and leading 
to the development of an occlusion. The fact that the introduction of the second- 
order closure boundary layer caused the seclusion to be nearly undetectable while 
the observations of maritime cold fronts clearly showed the presence of a seclusion 
suggests that there may be additional physical processes important in forming the 
seclusion that are not accounted for in the second-order closure simulation. 
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Figure A.l: Norwegiftn model schematic. 
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Figure A.2: SK model schematic. 


210 











Figure A.4: Croes section of potential temperature at hour 120 (AI). 
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Figure A.7: Cross section of potential temperatiire (K) at hour 108 (20C). 
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Figure A.8: 700 mb temper&ture distribution at hour 108 (20C). 
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